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Preface 


This  thesis  incorporates  the  various  methods  for  obtaining  the  unit 
risk  number  of  a  chemical  into  one  unified  process.  I  was  not  familiar 
with  the  chemical  risk  assessment  process  prior  to  working  on  this 
thesis.  My  interest  in  doing  research  on  the  chemical  risk  process  was 
mainly  founded  out  of  curiosity.  I  realized  the  decisions  associated 
with  hazardous  waste  are  enormous  and  if  I,  as  an  outsider,  could  shed  a 
little  light  on  just  a  small  piece  (unit  risk)  then  maybe  the  enormity 
of  the  task  could  be  reduced. 
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easy.  Joe  was  there  to  keep  me  in  line  and  always  had  positive  words  of 
encouragement.  Also,  I  thank  my  sponsors,  Lt.  Colonel  Clewell  and  Dr. 

Melvin  Andersen,  who  said  they  were  never  too  busy  for  ay  confused 
questions.  I  was  in  their  office  every  week  picking  their  brains,  and 
never  would  have  understood  the  chemical  risk  assessment  process  without 
their  help.  In  addition,  I  thank  my  reader,  Dr.  James  W.  Chrlssls,  for 
the  constructive  comments  which  have  made  this  thesis  more  palatable. 

And  finally,  I  thank  my  wife  Judy,  daughter  Theresa,  and  son  Paul,  who 
have  supported  me  through  this  entire  ordeal.  How  did  we  ever  find  the 
time  to  have  a  baby  son? 

Lastly,  I  wish  other  AFIT  students  will  be  as  fortunate  as  I. 

David  S.  Clement 
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Abstract 

For  a  given  chemical  there  are  usually  several  methods  for  estima¬ 
ting  the  risk.  Each  method  is  based  on  different  assumptions.  The 
arguments  are  plentiful  for  each  method,  but  which  method  best  estimates 
the  risk?  Choosing  one  method  over  another  could  lead  to  faulty  risk 
estimates,  thus  the  traditional  methods  have  been  very  conservative  to 
avoid  underestimating  the  risk.  With  advances  made  in  pharmacokinetics 
the  EPA  as  come  under  pressure  to  re-evaluate  its  procedure  for 
assigning  risk.  Which  method  or  methods  should  be  used  and  how  much 
emphasis  should  be  placed  on  each  one?  This  study  decomposes  the 
various  methods  into  their  corresponding  assumptions.  A  tree  diagram  is 
generated  to  describe  the  combinations  of  assumptions  that  make  up  each 
unit  rise  method,  a  subjective  weight  is  assigned  to  each  assumption 
(branch  of  the  tree)  to  characterize  its  validity  in  estimating  the 
risk.  From  this  a  weighted  average  of  risk  is  calculated.  A  procedure 
is  recommended  for  combining  expert  opinion  when  several  experts  are 
utilized  in  assigning  the  subjective  weights.  Two  examples  involving 
Methylene  Chloride  and  2,3,7,8-tetrachiorodibenzo-p-dioxln  illustrate 
the  decomposition  method  of  estimating  chemical  risk. 


DEALING  WITH  UNCERTAINTY  IN  CHEMICAL  RISK  ANALYSIS 


I.  Introduction 


Specific  Problem 

A  basic  problem  facing  the  U.S.  Environmental  Protection  Agency 
(EPA)  is  what  method  to  use  in  determining  cne  unit  risk  number  of  a 
chemical  that  have  shown  to  be  carcinogenic.  Traditionally, 
experimental  animal  studies  and  particular  statistical  methods  have  been 
used  in  determining  the  unit  risk  number.  Recently,  physiologically 
based  pharmacokinetic  models  (PBPK)  "have  been  developed  for  a  variety 
of  volatile  and  nonvolatile  chemicals,  and  their  ability  to  perform  the 
extrapolations  needed  in  risk  assessment  has  been  demonstrated”  (Clewell 
and  Andersen,  1985:111).  This  study  develops  an  analytical  method  that 
allows  for  several  models,  including  PBPK  models,  to  be  incorporated 
into  the  risk  assessment  process. 

Kev  Terms 

Carcinogen.  A  carcinogen  is  any  substance  that  is  associated  with 
producing  cancer: 

Cancer  is  now  considered  to  be  the  end  result  of  a  multistage 
process  in  which  a  large  number  of  endogenous  and  exogenous  factors 
interact,  simultaneously  or  in  sequence,  to  disrupt  normal  cell 
growth  and  division  (Wilkinson,  1987:844], 

Unit  Risk  Number.  "Risk  assessment  is  the  necessary  process  of 
determining  human  health  risks  from  exposure  to  chemicals"  (Menzel, 
1987:944).  The  unit  risk  number  of  a  chemical  is  the  probability  of  a 
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person  developing  cancer  if  exposed  to  a  set  concentration  cf  the 
chemical  in  the  air  breathed  or  in  the  food  and  water  that  is  ingested. 

A  chemical  may  have  several  unit  risk  numbers: 

1)  A  unit  risk  number  associated  with  continuous  inhalation  of 
a  specific  exposure  concentration  such  as  lug/m  . 

2)  A  unit  risk  number  associated  with  an  occupational  exposure 
such  as  being  exposed  to  1  part  per  million  (ppm),  8  hours  per  day, 
5  days  per  week,  and  48  weeks  per  year. 

3)  Or  a  unit  risk  number  associated  with  ingesting  the 
chemical,  such  as  when  drinking  contaminated  ground-water.  The 
unit  risk  number  would  then  be  a  risk  per  1  mg/liter. 

For  example:  suppose  the  unit  risk  number  of  methylene  chloride 

(DCM),  an  important  solvent  used  in  paint  removers,  is  4.7xl0-7  for 

continuous  inhalation  of  1  microgram  of  DCH  per  cubic  meter  of  air  per 

day  (1  ug/m3,  or  lxlO6  g/m3).  Then  a  person  breathing  1  ug  of 

methylene  chloride  per  cubic  meter  of  air  each  day  over  his  life  time 

would  have  a  4.7x10"'  chance  of  developing  cancer.  Likewise,  the  higher 

the  concentration  the  higher  the  probability.  The  unit  risk  number  is 

usually  given  as  the  probability  of  developing  cancer  from  continuous 

exposure  to  1  ppm,  1  mg/kg/day,  or  1  pg/kg/day  of  chemical. 

Pharmacokinetics.  Pharmacokinetics  is  "the  quantitative  study  of 

the  metabolic  processes  of  absorption,  distribution,  biotransformation, 

and  elimination"  (Calabrese,  1987:619).  Physiologically-based 

pharmacokinetic  models  (PBPK)  ere  defined  as  follows: 

They  are  essentially  mechanistic  models  that  try  to  account 
quantitatively  over  time  for  the  various  pharmacokinetic  processes 
that  involve  an  agent  [chemical]  of  concern  from  the  time  the  agent 
reaches  a  site  of  absorption  to  the  time  an  interaction  occurs 
between  the  agent,  its  metabolites,  and  various  body  tissues 
[Calabrese,  1987:620]. 

PBPK  models  are  designed  to  determine  the  relationship  between  the 
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internal  dose  in  the  tissues  and  the  amount  of  applied  dose  (EPA,  July 
1987:47). 

General  Background 

The  Risk  Assessment  Process.  The  EPA  is  responsible  for  assigning 

unit  risk  numbers  to  chemicals  that  have  been  shown  to  be  carcinogenic. 

A  chemical  is  considered  carcinogenic  if  it  has  been  found  to  promote 

cancer  in  animals  during  long-term  experiments  (EPA,  September 

1986:33995).  Usually,  the  EPA  contracts  with  a  lab  to  test  a  chemical 

that  is  suspected  to  be  carcinogenic  as  defined  by  EPA  Guidelines  for 

Carcinogen  Risk  Assessment  (EPA,  September  1986).  The  EPA  then 

evaluates  the  laboratory  studies  ... 

according  to  sound  biological  and  statistical  considerations  and 
procedures.  Results  and  conclusions  concerning  the  agent 
[chemical]  derived  from  different  types  of  information  ...  are 
melded  together  into  a  weight-of-evidence  determination.  The 
strength  of  evidence  supporting  a  potential  human  carcinogenicity 
judgment  is  developed  in  a  weight-of-evidence  stratification  scheme 
[EPA,  September  1986:33994]  ... 

which  is  used  in  assigning  a  unit  risk  number  to  the  chemical.  The  EPA 
uses  more  than  the  laboratory  results  in  determining  the  unit  risk 
number.  The  EPA  asks  for  comments  from  the  scientific  community.  Both 
industrial  and  environmental  groups  present  arguments  for  lowering  or 
raising  the  unit  risk  number.  Based  on  the  experimental  findings  and 
community  debate,  the  EPA  then  publishes  a  unit  risk  number  for  the 
chemical . 

The  Armstrong  Aerospace  Medical  Research  Laboratory,  Toxic  Hazards 
Division,  is  responsible  for  advising  on  the  toxicity  of  AF  chemicals 
and  has  provided  input  to  the  EPA  on  issue?  relating  to  AF  use  of  toxic 
chemicals.  The  AF  is  tasked  with  cleaning  up  their  hazardous  waste 
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sites,  and  spends  approximately  S00  Million  dollars  a  year  on  cleaning 
up  their  hazardous  waste  sites.  The  unit  risk  number  dictates  the 
cleanliness  standard  the  AF  must  aeet.  A  high  probability  unit  risk 
nuaber  neans  tougher  exposure  standards,  and  tougher  standards  Bake  for 
higher  clean  up  costs,  and  the  higher  costs  add  up  to  fewer  hazardous 
waste  sites  getting  cleaned.  With  increasing  public  support  for 
cleaning  all  hazardous  waste  sites,  the  AF  is  driven  to  use  its  limited 
budget  for  cleaning  up  as  many  hazardous  waste  sites  as  possible.  If 
the  unit  risk  number  is  unreasonably  conservative  the  budget  is  wasted 
on  over  cleaning.  But  if  the  unit  risk  is  set  too  low,  people  are 
exposed  to  a  greater  risk  than  is  publicly  acceptable.  Again  society 
bears  the  cost  of  the  medical  expenses  in  treating  the  exposed 
individuals.  Of  course  nobody  wants  an  unsafe  environment,  but  wasting 
precious  resources  (budget)  on  over  stringent  clean  up  standards  is  poor 
policy.  Good  science  should  lead  to  sound  policy  that  all  of  society 
can  live  with. 

Statistical  Methods.  The  traditional  method  of  assessing  risk,  as 
stated  earlier,  starts  with  long-term  animal  exposure  studies.  During 
these  tests,  animals  (usually  mice)  are  exposed  to  various  high 
concentrations  of  the  chemical  in  their  air,  or  water,  or  feed.  If  the 
animals  develop  tumors,  then  the  chemical  is  classified  as  a  probable 
human  carcinogen.  The  probability  of  the  animals  developing  cancer  from 
various  exposure  doses  is  determined  and  is  used  for  extrapolating  the 
probability  of  developing  cancer  in  humans. 

The  statistical  methods  for  determining  unit  risk  numbers  have  been 
Improved  over  the  years,  but  they  still  rely  on  extrapolating  the  risk 
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in  animals  to  the  risk  in  humans.  The  methods  have  assumed  humans  to  be 
more  sensitive  to  the  tested  chemical  than  the  lab  animals  (body  surface 
method).  This  may  not  always  be  true  (Calabrese,  1987:620).  Some 
methods,  when  extrapolating  from  high  dose  to  low  dose  (dose  response 
models),  ignore  the  possibility  that  below  a  certain  threshold  level  the 
chemical  may  pose  no  risk  of  cancer  at  all.  In  addition,  bioassays  are 
performed  exposing  the  animal  to  contaminated  air  and  water,  for  which 
the  concentrations  are  known  (external  factors),  and  the  internal  tissue 
concentrations  of  the  chemical  are  based  on  these  external 
concentrations  and  the  animal's  breathing  rate.  The  chemical  is  assumed 
to  be  completely  absorbed  into  the  animal's  body,  which  completely 
neglects  the  physiological  aspects  of  the  body  in  dealing  with  the 
chemical  ( Henze 1,  1987:944).  Every  year  more  and  more  research  is  done 
that  adds  to  our  understanding  of  human  and  animal  physiology.  The  EPA 
wants  to  take  advantage  of  this  information  when  assigning  unit  risk 
numbers. 

PBPK  Methods.  During  the  past  few  years,  there  have  been  great 
strides  made  in  understanding  the  physiological  basis  of  how  mammals 
process  chemicals  in  the  blood  stream. 

By  using  physiologically  based  mathematical  models  of  the 
transport,  distribution,  and  metabolism  of  toxic  chemicals, 
scientists  can  predict  and  confirm  the  organ  or  intracellular  dose 
of  the  chemical  by  experimental  analysis  in  exposed  animals.  By 
this  process  the  exposure  concentration  is  converted  to  a  tissue 
dose  [Menzel,  1987:945]. 

As  with  most  new  methods,  these  PBPK  approaches  have  yet  to  gain 
widespread  acceptance.  The  advantage  of  PBPK  models  over  existing 
statistical  methods  is  that  PBPK  models  are  "based  to  a  large  extent  on 
the  actual  physiology  of  the  organism  [Clewell  and  Andersen,  1965:114]." 
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"Since  the  relationship  between  the  effective  dose  [internal  dose, 
target  dose,  or  surrogate  dose]  and  the  administered  dose  [external 
dose]  m ay  be  nonlinear,  pharmacokinetic  considerations  can  significantly 
modify  low-dose  risk  estimates  regardless  of  the  model  employed  [Federal 
Register,  1985:81]."  PBPK  models  have  provided  a  better  understanding 
of  the  way  the  body  processes  chemicals,  but  as  yet  have  only  been 
marginally  used  in  risk  assessment  (Andersen  and  Clewell,  1988). 

The  EPA  is  interested  in  using  new  approaches  for  assigning  the 
unit  risk  number  (Shabecoff,  1968).  The  PBPK  models  can  provide  more 
realistic  data  on  the  effects  chemicals  have  on  the  body.  Four 
alternatives  have  been  suggested  where  PBPK  models  may  aid  the  EPA  in 
assigning  the  unit  risk  number  (Cohn,  1987:2): 

1)  PBPK  models  are  ignored  (too  new  and  not  proven). 

2)  PBPK  models  could  assist  in  extrapolating  from  the 
high  to  low  dose. 

3)  They  could  assist  in  the  species-to-species  extrapolation. 

4)  PBPK  models  could  be  used  in  both  dose  and  species 
extrapolations. 

Problem  Statement 

The  EPA  is  interested  in  improving  the  way  the  unit  risk  numbers 
are  determined  for  carcinogens.  What  is  needed  is  a  method  that 
addresses  the  uncertainties  associated  with  the  various  techniques  for 
assigning  chemical  risk  and  to  combine  them,  including  the  PBPK  modeling 
technique,  into  one  cohesive  unit  risk  assessment  process. 

Scone.  This  study  will  develop  a  methodology  for  assigning  unit 
risk  numbers  by  combining  both  the  traditional  methods  and  the  PBPK 
modeling  techniques  into  one  risk  assessment  process,  with  emphasis  on 
methodology  for  assigning  the  unit  risk  number  and  not  on  the  unit  risk 
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number  itself.  The  study  is  limited  to  two  chemicals:  methylene 
chloride  (DCM)  and  2,3,7,8-tetrachlorodibenzo-p-dioxin  (dioxin).  Both 
chemicals  are  used  for  illustrating  the  method. 

Methodology.  Several  methods  can  be  evaluated  for  determining  the 
risk  of  a  chemical  by  decomposing  them  into  their  corresponding 
assumptions.  A  subjective  probability  is  then  assigned  to  each 
assumption  to  characterize  its  validity  in  estimating  the  risk.  From 
this  a  weighted  average  of  risk  is  calculated. 

A  general  discussion  of  dealing  with  uncertainty  of  the  assumptions 
using  tree  diagrams  is  presented  first.  Then  how  to  combine  experts' 
subjective  weights  is  reviewed.  And  finally  two  examples  involving  DCM, 
an  important  industrial  solvent  used  by  the  AF  for  paint  striping 
operations,  and  dioxin  will  illustrate  the  process.  DCM  shows  how  to 
incorporate  pharmacokinetics  into  the  risk  assessment  process  since  it 
is  one  of  the  few  chemicals  to  have  a  good  PBPK  model  built  for  it.  And 
dioxin  was  chosen  to  show  the  flexibility  of  the  decomposition  process 
because  the  risk  associated  with  it  was  developed  differently  than  DCM. 
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II.  Model  Development 


There  are  a  variety  of  ways  for  determining  the  risk  of  a 
carcinogen.  Chapter  I  discussed  the  overall  process.  Here  the  process 
will  be  examined  in  a  flowchart  like  fashion.  Hogarth  (1986:146) 
suggests  "decomposing  the  problem  into  appropriate  components",  thus  it 
becomes  easier  to  understand. 

Flowchart 

The  unit  risk  number  assessment  flowchart  is  shown  in  Figure  1. 

Each  sequence  in  the  flowchart  is  as  follows: 

Bioassav.  The  first  step  in  assessing  the  unit  risk  number  for  any 
suspected  carcinogen  is  to  conduct  a  bioassay  to  determine  the  effect  it 
has  upon  animals.  The  bioassays  are  laboratory  experiments  carried  out 
on  animals,  usually  rats  or  mice,  in  the  expectation  of  determining  the 
risk  (probability  of  developing  cancer,  or  some  other  measured  effect) 
associated  with  exposure  to  the  suspected  carcinogen.  The  risk 
associated  with  the  bioassay  is  dependent  on  several  factors: 

1)  Suspected  Carcinogen,  (in  this  case  dioxin  and  DCH). 

2)  Dose,  the  concentration  administered  to  the  animal. 

3)  Route,  exposure  via  air,  water,  food,  or  touch. 

4)  Effect,  does  the  suspected  carcinogen  cause  cancer, 
tumors,  disrupt  the  immune  system  or  reproductive  track,  alter  cell 
growth,  or  maybe  create  nervous  system  disorders? 

5)  Species,  several  different  species  could  be  used,  but  some 
animals  do  not  react  the  same  way  as  a  human  would  if  exposed  to 
the  same  chemical.  Depending  on  what  effect  is  expected  a  specific 
animal  is  used. 

The  administered  dose  is  usually  very  high  in  order  to  promote  the 
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Figure  1.  Chemical  Risk  Assessment  Flowchart. 


development  of  the  effect  In  the  animal.  If  the  dose  is  too  small,  the 
effect  may  not  develop  prior  to  the  animal  dying  a  natural  death.  Once 
the  bioassay  is  complete,  it  can  be  used  in  several  ways  to  extrapolate 
the  risk  from  animal  to  human. 

Species  Internal  Dose.  The  species'  internal  dose  is  calculated 
in  one  of  two  ways:  using  the  applied  dose  technique  or  the  PBPK  model. 

1)  Applied  dose  technique.  The  animal’s  administered 

dose  (ppm)  is  converted  to  an  Internal  dose.  The  administered  dose 
is  multiplied  by  the  animal's  breathing  rate  and  divided  by  its 
body  weight  to  obtain  an  internal  dose  in  units  of  mg/kg/day.  For 
example,  if  a  mouse  is  exposed  to  4000  ppm  of  DCM  for  6  hours  a 
day,  5  days  a  week,  then  the  average  daily  internal  dose  is: 

int.  dose  *  (4000  ppm)  *  (3.4766  mg/m3/l  ppm) 

*  (0.043  m3/day>  *  (1/0.0345  kg) 

*  (6  hrs/24  hrs)  *  (5  days/7  days) 
int.  dose  *  3095  mg/kg/day 

2)  PBPK  model  technique.  This  time  the  administered  dose  is 
converted  to  an  internal  dose  (tissue  dose  or  dose  surrogate)  using 
a  physiologically-based  pharmacokinetic  (PBPK)  model  (Andersen  et 
al,  1987).  Refer  to  appendix  A  for  details  on  the  PBPK  model. 

Risk-per-Dose.  The  risk-per-dose  is  obtained  by  taking  the 

bioassay  results  and  the  species  internal  dose  and  plugging  them  into  a 

dose  response  model.  One  of  these  dose  response  models  is  the  linear 

multistage  (LMS)  model  which  provides  a  95  percent  lower  confidence 

limit  (LCL)  on  dose.  The  LMS  model  has  been  commonly  used  by  the  EPA 

since  it  provides  a  conservative  upper  confidence  limit  on  risk  (Munro 

and  Krewakl,  1981:553;  EPA,  November  1987:27-31).  Another  method  uses 

the  maximum  likelihood  estimator  (MLE)  for  dose  instead  of  the  95 

percent  LCL.  Traditionally,  the  dose  response  risk  has  been  obtained 

using  safety  factors.  The  risk  associated  with  the  bioassay  is 
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multiplied  by  some  safety  factor  (e.g.,  100  -  1000)  to  obtain  the  risk 

for  a  given  dose  (EPA,  November  1987:34).  For  example,  using  the  LMS 

model  the  number  of  animals  which  developed  cancer  for  the  given  sample 

size  and  exposure  rate  is  extrapolated  to  a  95  percent  lower  confidence 

-8 

limit  dose  for  a  1x10  chance  of  cancer.  Refer  to  appendix  B  for 
details  on  the  LMS  dose  response  model. 

Human  Internal  Dose.  The  human  internal  dose  is  calculated  much 
the  same  way  as  the  species  internal  dose,  but  this  time  the  human 
exposure  dose  is  set  to  some  standard  setting,  such  as  a  continuous 
exposure  to  1  ug/m 3.  The  human  internal  dose  is  calculated  in  one  of 
two  ways:  using  the  applied  dose  technique  or  the  PBPK  model. 

1)  Applied  dose  technique.  The  exposure  of  interest  is  then 
converted  to  an  internal  human  dose  as  follows: 

int.  dose  »  (1  ug/m3)  *  (.001mg/l  ug)  *  (20  m3/day)  *  (1/70  kg) 
int.  dose  *  0.000286  mg/kg/day 

2)  PBPK  model  technique.  In  this  case  the  dose  of  interest 
is  converted  using  the  physiologically-based  pharmacokinetic  (PBPK) 
model  (Andersen  et  al,  1987). 

Target  Dose.  The  human  internal  dose  is  multiplied  by  the  species- 
to-human  pharmacodynamics  correction  factor  to  get  the  tissue  target 
dose.  There  are  two  correction  factors: 

1)  Body  weight  (8W).  Since  the  internal  dose  is  already  in 
units  of  weight  (kg)  the  BW  correction  factor  is  1. 

2)  Body  surface  (BS).  When  using  the  body  surface  factor  a 
human  is  considered  more  sensitive  to  the  chemical  exposure  than 
the  animal.  The  internal  dose  calculated  above  is  multiplied  by 
the  body  surface  ratio  to  obtain  a  larger  internal  dose.  The  body 
surface  ratio  is  the  human's  weight  to  animal's  weight  raised  to 
the  l/3rd  power.  For  a  mouse  the  BS  correction  factor  is: 

BS  =  (70  kg/0.0345  kg)1/3 

BS  *  12.66 
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Thus,  the  body  surface  internal  dose  is  approximately  12.7  times 
larger  and  will  make  the  unit  risk  number  of  a  chemical  12.7  times 
higher  than  if  calculated  using  the  body  weight  method. 

Unit  Risk  Humber.  After  obtaining  the  animal's  risk-per-dose 

factor  and  human  target  dose,  the  unit  risk  number  for  a  given  exposure 

can  be  determined.  For  example: 

Unit  Risk  =  risk-per-dose  *  target  dose 

-8 


Unit  Risk 


1x10 


1.3E-5  ag/kg/day 


*  2.8xl0-4  mg/kg/day 


.-7 


Unit  Risk  =  2.15x10  for  a  continuous  exposure 
to  1  ug/kg/day  of  chemical 

The  unit  risk  number  could  either  be  the  probability  of  developing  the 

effect  at  some  set  dose,  as  in  the  above  example,  or  it  could  be  the 

maximum  recommended  dose  (risk-specific  dose)  at  a  set  probability 
-4  -5  -6 

(e.g.,  10  ,10  ,  or  10  ),  as  is  the  case  for  dioxin. 


Tree  Diagram 

There  has  been  much  controversy  on  the  appropriate  method  for 
determining  the  unit  risk  number  of  a  chemical.  Any  policy  put  forth 
that  suggests  a  specific  unit  risk  number  for  a  particular  chemical  is 
usually  criticized  for  some  of  its  assumptions.  There  is  always  someone 
who  disagrees  on  how  the  assumptions  were  made  or  what  techniques  were 
utilized.  In  a  majority  of  cases  the  unit  risk  number  is  developed 
using  a  combination  of  the  most  likely,  and  sometimes  most  conservative, 
assumptions  for  a  given  chemical  and  effect.  The  problem  with  this 
approach  is  that  other  assumptions,  even  though  less  likely,  are  not 
considered  in  the  final  unit  risk. 

If  all  potential  assumptions  could  be  weighted  according  to  their 
likelihood  of  estimating  the  unit  risk  number,  then  a  weighted  average 
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of  the  unit  risk  numbers  could  be  determined.  In  other  words,  by  taking 
the  various  methods  for  determining  the  unit  risk  number  and  weighting 
each  according  to  the  assumptions  used,  then  a  weighted  average  unit 
risk  number  can  be  calculated,  where  the  weights  reflect  the  decision 
makers  belief  in  each  assumption's  ability  to  accurately  determine  the 
unit  risk.  The  assignment  of  the  unit  risk  number  could  then  take  on 
the  form  of  a  tree  diagram  with  the  branches  representing  each  potential 
assumption;  at  the  end  of  the  branches  would  be  the  unit  risk  numbers 
for  the  various  combinations  of  assumptions.  Weights  would  be  assigned 
by  the  decision  maker  (EPA)  to  each  branch  for  the  likelihood  of  that 
branch  in  estimating  the  unit  risk  number.  Thus,  all  arguable 
combinations  of  assumptions  could  be  considered  in  obtaining  the  unit 
risk. 

The  weights  fall  into  two  classes:  preference  weights  and 
likelihood  weights.  The  weights  associated  w'th  the  species 
pharmacokinetics,  human  pharmacokinetics,  and  pharmacodynamics  reflect 
the  decision  makers  belief  (preference)  of  how  well  the  branch  predicts 
the  unit  risk  of  a  chemical.  The  values  of  the  PBPK  variables  and  their 
respective  probabilities  (likelihood  weights)  reflect  the  likelihood  of 
the  actual  or  true  value  of  the  variable  (a  percentile  of  the  variable's 
probability  distribution  on  value).  Nixing  the  preference  weights  with 
likelihood  weights  has  no  theoretical  basis,  so  the  decomposition  method 
splits  the  tree  diagram  into  two  areas:  preference  weights  and 
likelihood  weights  (probabilities).  The  likelihood  weights  associated 
with  the  PBPK  variables  are  grouped  into  one  tree  and  an  expected  value 
of  unit  risk  is  calculated.  Then,  the  expected  value  of  unit  risk  is 
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Inserted  In  the  preference  weight  tree  and  a  weighted  average  of  unit 
risk  is  determined. 

There  are  two  major  drawbacks  in  using  every  tree  combination  for 
assigning  the  unit  risk  number:  1)  there  is  not  enough  time  or  money  to 
study  every  potential  combination,  and  2)  how  to  get  a  consensus  on  what 
weight  to  assign  each  branch.  Obtaining  consensus  for  the  weights  will 
be  discussed  in  the  next  chapter,  and  by  only  developing  a  tree  based  on 
current  information  will  the  need  to  study  every  potential  combination 
be  eliminated.  As  more  information  is  learned  on  a  specific  chemical 
the  tree  may  have  to  be  enlarged  to  accommodate  this  new  information. 

Hew  weights  would  be  assigned  to  both  new  and  existing  branches  and 
another  weighted  average  unit  risk  number  calculated.  Hence,  the  tree 
diagram  would  be  the  framework  for  analyzing  the  unit  risk  of  a 
chemical.  Dioxin  and  DCH  are  used  to  illustrate  building  trees  for 
calculating  chemical  risk. 

Dioxin  Tree.  The  information  used  for  developing  a  dioxin  tree 
came  from  "A  Cancer  Risk-Specific  Dose  Estimate  For  2,3,7,8-TCDD”  from 
the  EFA  (EPA,  November  1987).  This  report  contains  a  summary  of 
acceptable  studies  (14  in  all)  for  generating  a  risk-specific  dose  for 
dioxin,  which  is  the  acceptable  dose  for  a  unit  risk  number  of  1x10  6. 
The  studies  were  grouped  according  to  assumption  in  developing  a  tree 
diagram.  Four  studies  used  the  traditional  safety  factor  approach, 
sometimes  referred  to  as  the  no  effect  level  (HOEL)  approach,  meaning 
that  below  some  threshold  value  the  risk  is  assumed  to  be  zero.  The 
remaining  ten  studies  assumed  no  threshold  and  thus  there  is  a  chance  of 
developing  the  effect,  in  this  case  cancer,  for  even  the  smallest  doses. 
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These  ten  studies  were  further  categorized  according  to  the  dose 
response  model  used  in  extrapolating  the  dose  for  a  given  risk.  The 
tree  for  dioxin  is  shown  in  Figure  2. 

PCM  Tree.  The  tree  diagram  generated  for  DCM  is  more  symmetrical 
than  the  dioxin  tree  (Figures  3  and  4).  Most  of  the  information  for 
developing  the  DCM  tree  diagram  came  from  Andersen  and  Clewell  (1988). 
Figure  3  contains  the  major  assumptions  employed  by  federal  agencies  in 
determining  a  unit  risk  number.  Figure  4  further  defines  the  PBPK 
branch  of  Figure  3  by  taking  the  mouse  and  human  pharmacokinetic  models 
and  dividing  them  into  several  PBPK  variables.  According  to  the  EPA 
(EPA,  July  1987:73)  the  PBPK  predictive  model  of  Clewell  and  Andersen 
for  OCM  did  not  take  into  consideration  several  alternative  ways  of 
predicting  DCM's  unit  risk  number,  such  as  the  chemical  using  a 
different  pathway  to  affect  the  body's  tissues.  Also  the  PBPK  model  did 
not  allow  for  uncertainty  of  its  parameter  values.  This  tree  allows 
those  alternatives  to  be  included  in  the  assessment  process  and  weighted 
as  to  their  ability  to  predict  the  unit  risk  for  DCM.  Thus,  by  varying 
the  parameter  values,  the  uncertainty  in  which  value  to  use  can  be 
addressed,  and  a  much  broader  and  more  believable  approach  generated  for 
determining  DCM's  unit  risk  number. 

The  DCM  tree,  Figure  3,  is  divided  into  three  sections: 

1)  Dose  response  model,  which  is  modeled  by  the  LMS  model. 

No  other  model  has  been  used  to  obtain  the  unit  risk  number  for 
DCM. 

2)  Species-to-human,  which  is  broken  into  two  subsections: 

A)  Species  pharmacokinetics,  which  was  modeled 
traditionally  using  the  applied  dose  method  and  only  during 
the  last  couple  of  years  has  the  PBPK  model  been  considered. 
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Figure  3.  Methylene  Chloride  (DCM) 


Tree  Diagram. 
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B)  Pharmacodynamics,  which  is  modeled  using  body  surface 
and  body  weight  calculations. 

3)  Hunan  pharmacokinetics,  which,  as  with  section  2A  above, 
has  been  traditionally  modeled  using  the  applied  dose  method  and 
only  during  the  last  several  years  has  been  evaluated  using  the 
PBPK  model. 

PCM  PBPK  Tree.  The  PBPK  portion  of  the  DCM  tree,  Figure  4,  is 
broken  into  two  sections: 

1)  Pathway,  the  metabolism  occurs  via  two  pathways:  "one 
dependent  on  oxidation  by  mixed  function  oxidases  (HFO)  and  the 
other  dependent  on  glutathione  S-transferases  (GST)  [Andersen  et 
al.,  1987}",  or  DCM  could  act  directly  on  the  tissues  of  concern. 
Thus  the  tree  is  broken  into  three  additional  assumptions  regarding 
the  chemical  mechanism  responsible  for  causing  cancer: 

A)  GST  pathway, 

B)  MFO  pathway, 

C)  DCM  directly. 

2)  Biochemical  constants  (metabolic  constants  and  kinetic 
constants)  which  are  further  broken  into: 

A)  Al  and  A2,  relative  activity  ratios  of  lung  enzymes 
to  liver  enzymes  in  the  MFO  and  GST  pathways  respectively. 

B)  K.,  first-order  rate  constant  for  metabolism  of  DCM 
by  the  GSTpathway  in  the  liver. 

C)  ,  the  maximum  velocity  of  metabolism  by  the  MFO 
pathway  in  cfie  liver. 

The  choice  of  limiting  the  biochemical  constants  to  the  four  listed 
above  (Al,  A2,  K.,  V  )  was  made  by  Dr.  Andersen. 

I  Max 

Assumptions.  The  National  Toxicology  Program’s  (NTP) 
Inhalation  bloassay  (NTP,  1986)  found  that  DCM  caused  cancer  in  male  and 
female  mice.  With  only  one  bloassay  used  for  analysis  there  is  no  need 
for  a  bioassay  branch  within  the  DCM  tree.  Also  there  is  not  a  branch 
for  the  dose  response  model  node,  since  only  the  LMS  dose  response  model 
has  been  used  by  federal  agencies  in  determining  DCM’s  unit  risk  number. 
All  other  PBPK  values,  other  than  Al,  A2,  K.,  and  V  ,  are  set  to 

C  BflX 
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nominal  values. 


For  a  discussion  of  these  values  see  appendix  A. 


Model  Summary 

The  various  ways  of  obtaining  the  risk  for  dioxin  and  OCM  have  been 
decomposed  into  their  different  assumptions.  The  tree  diagram  allows 
the  decision  maker  to  focus  on  each  assumption  and  assign  a  weight  as  to 
his  belief  that  the  assumption  will  aid  in  providing  a  realistic  risk. 
The  decision  maker  will  most  likely  seek  expert  advice  in  determining 
the  weights  of  each  branch,  and  would  most  likely  consult  with  more  than 
one  person.  Combining  the  experts'  weights  into  one  weight  for  each 
branch  poses  several  problems  which  will  be  discussed  in  the  next 
chapter. 
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III.  Combining  Expert  Opinion 


The  Problem 

The  tree  diagrams  developed  in  the  previous  chapter  required 
preference  weights  to  be  assigned  to  the  various  branches.  A  decision 
maker  nay  not  know  all  the  facts  about  a  certain  chemical  and  the 
different  models,  assumptions,  and  techniques  that  have  been  developed 
for  predicting  the  unit  risk.  So  a  decision  maker  may  seek  guidance, 
that  is,  an  opinion  from  an  expert,  someone  who  is  familiar  with  the 
chemical  in  question,  in  assigning  these  weights.  In  the  case  of 
seeking  only  one  expert's  opinion  about  a  branch's  weight,  as  in  this 
thesis,  the  expert's  weighting  becomes  the  decision  maker's  weight.  A 
problem  arises  when  the  decision  maker  asks  more  than  one  expert  for  an 
opinion  in  assigning  preference  weights  to  the  branches.  Different 
experts  may  have  different  weights,  and  combining  these  weights  is  not  a 
simple  matter. 

Combining  Techniques 

An  expert's  preference  weight  for  each  branch  can  be  thought  of  as 
his  subjective  probability;  so,  combining  a  group  of  experts'  weights  is 
like  combining  subjective  probabilities.  There  are  two  schools  of 
thought  for  combining  subjective  probability  (Seaver,  1976:27;  Seaver, 
1978:3;  Edwards  et  al.,  1979:2): 

1)  Mathematical  aggregation  models  -  which  use  some 
mathematical  technique  for  averaging  each  expert’s  opinion. 

2)  Group  behavioral  techniques  -  which  uses  interaction 
procedures  for  developing  a  group  consensus  from  the  individuals' 
opinions. 
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Mathematical  Approach.  There  are  numerous  methods  for 
mathematically  combining  subjective  probability,  which  include  the 
arithmetic  average,  the  geometric  average,  the  ranked  or  weighted 
average,  the  bayeslan  and  pseudo  bayesian,  and  the  personal  interview. 
This  is  not  all  inclusive  of  the  methods  available  for  combining  expert 
opinion. 

Arithmetic  Average.  The  arithmetic  average  is  defined  as  the 
average  of  the  different  expert's  weights: 

v  1 

x  =  g  *  SUM(Xl )  for  i  *  1  to  N  (1) 

where,  Xj  =  the  ith  expert's  weight 

and,  H  *  the  number  of  experts 

The  arithmetic  average  may  also  take  on  the  trimmed  form  where  the  high 
and  low  expert  weights  are  discarded.  Depending  on  the  number  of 
experts  the  two  highest  and  lowest  expert  weights  may  be  discarded.  In 
this  way  the  weights  only  include  the  moderate  values  and  any  extremes 
are  eliminated.  The  trimmed  arithmetic  average  is  defined: 

*t  3  fib  *  SUM<V  for  i  *  1  to  N-J  (2) 

where,  j  *  2,  for  the  one-trimmed  case  when  the  high  and  low 
value  are  discarded. 

or,  j  *  4,  for  the  two-trimmed  case  when  the  two  highest 
and  two  lowest  values  are  discarded. 

Geometric  Average.  The  geometric  average  is  similar  to  the 

mathematical  average,  but  the  logarithm  of  the  experts'  weights  are 

averaged: 

Log(X)  =  J  *  SUM(Log(Xi ) )  for  i  *  1  to  H  (3) 

Ranked  or  Weighted  Average.  The  ranked  (weighted)  average  is 
used  when  the  decision  maker  does  not  believe  each  expert's  weight 
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should  be  given  equal  treatment.  The  decision  maker  may  have  some 

reason  for  believing  that  one  expert's  opinion  may  be  more  credible  than 

that  of  another.  In  this  case  the  decision  maker  may  rank  the  experts 

from  least  credible  to  most  credible  or  the  decision  maker  may  ask  the 

experts  to  rank  themselves  or  one  another.  The  ranked  (weighted) 

average  is  defined  as  follows: 

*r  =  SUM(WjXj )  for  i  =  1  to  N  (4) 

where,  Wi  =  I^/SUMO^) 

and,  Rj  *  rank  of  the  ith  expert 

Bayesian.  The  Bayesian  approach  is  considered  the  most 

realistic  way  to  combine  expert  subjective  probabilities  because  a 

decision  maker's  weight  for  a  branch  is  conditioned  upon  his  prior 

Information  and  upon  each  expert's  opinion.  Also,  each  expert's  opinion 

is  conditioned  upon  his  prior  information  (Morris,  1974:1235;  Hogarth, 

1986:145).  Quantifying  these  prior  probabilities  is  extremely  difficult 

in  all  but  the  simplest  cases  (Morris,  1974:1236;  Hogarth,  1975:283). 

Pseudo-Bayesian.  Since  the  Bayesian  approach  is  difficult  to 

apply,  an  alternative  has  been  proposed: 

If  one  expert's  distribution  of  the  unknown  state  of  nature  is 
treated  as  a  prior  and  the  second  expert's  distribution  as  a 
likelihood,  then  Bayes'  theorem  provides  a  mechanism  for  combining 
them.  This  procedure  can  then  be  repeated  successively  to  combine 
the  opinions  of  all  the  experts  (Chatter jee,  1987:278]. 

With  this  approach  the  decision  maker  does  not  have  to  assess  each 

expert's  prior  information. 

Interaction  Approach.  There  are  alternatives  for  combining  expert 
opinion.  These  methods  are  not  mathematically  oriented,  in  fact  they  do 
not  combine  each  expert's  opinion,  but  instead  use  interaction 
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techniques  so  as  to  Induce  an  atmosphere  for  free  exchange  of 
information  in  the  hope  a  consensus  of  opinion  may  develop.  The  Delphi 
method  is  a  popularly  accepted  method  for  seeking  opinions  from  a  group 
of  experts. 

Delphi  Method.  ( Chatter jee,  1987:285-286)  The  Delphi  method 
uses  an  iterative  approach  for  developing  a  consensus  of  opinion.  A 
survey  is  conducted  where  a  questionnaire  is  distributed  to  the 
participants  (experts).  The  experts  are  asked  not  to  personally  contact 
one  another  as  to  why  they  have  a  particular  opinion.  After  the 
questions  are  completed  the  questionnaires  are  collected  and 
statistically  compiled  and  the  results  are  made  known  to  the 
participants.  Who  are  again  asked  to  answer  the  questions  given  the 
results  from  the  first  iteration,  but  during  the  second  iteration  the 
participants  are  asked  to  justify  their  answers.  During  the  third 
iteration  the  participants  are  given  each  others  reasons  for  their 
opinions.  The  questions,  justification,  and  results  are  repeated  a 
fourth  time.  This  iterative  procedure  allows  each  participant  to  re¬ 
evaluate  his  opinion  three  or  four  times  with  the  others  participants' 
reasons  and  opinions.  The  Delphi  method  may  lead  to  a  consensus  of 
opinion,  but  a  consensus  is  not  forced. 

Applying  the  Delphi  method  to  the  tree  weights  would  allow  each 
expert  to  assign  an  opinion  (weight)  to  each  branch  of  the  tree.  During 
the  first  iteration  the  experts  would  only  have  their  own  prior 
knowledge  and  experience  to  call  on.  The  first  iteration  weights  would 
then  be  made  available  to  all  the  experts.  The  experts  would  then  be 
given  a  chance  to  re-evaluate  their  weights  and  to  state  their  reasons. 
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The  second  iteration's  results  with  reasons  are  Bade  available  and  a 
third  and  fourth  iteration  may  be  performed.  In  this  way  each  expert  is 
permitted  to  update  his  weight  given  the  new  prior  information. 

On  the  positive  side,  the  Delphi  method  (Dalkey,  1967:3): 

1)  provides  anonymity  and  thus  keeps  a  individual  from 
intimidating  others. 

2)  provides  controlled  feedback  which  can  be  used  for 
updating  the  expert's  opinion. 

3)  also  provides  a  statistical  group  response  so  that  each 
expert  knows  where  he  stands. 

On  the  negative  side,  the  responses  and  reasons  generated  by  each  expert 
could  become  less  convincing  with  each  iteration  since  the  expert  may 
grow  tired  of  the  procedure. 

Face-to-Face.  The  Delphi  method  does  not  allow  for  personal 
Interaction  and  is  limited  as  to  how  thorough  the  feedback  between 
experts  is  with  each  iteration.  If  the  experts  could  openly  confront 
and  discuss  their  beliefs  and  understandings,  then  a  consensus  of 
opinion,  if  possible,  may  more  easily  be  formed.  A  drawback  to  the 
face-to-face  interaction  is  that  an  individual  could  dominate  or  be 
intimidated  by  others,  thus  forcing  others  to  accept  or  reject  his 
opinion.  A  consensus  of  opinion  is  ideal,  but  only  if  it  is  not  forced. 
Another  problem  is  that  some  people  are  just  stubborn  because  they  don't 
want  to  be  thought  of  as  the  one  who  compromised.  This  is  sometimes 
referred  to  as  hardening  (Chatter jee,  1987:286-287).  Thus,  for  the  sake 
of  pride,  a  consensus,  if  possible,  falls  to  be  reached.  On  the 
positive  side,  by  bringing  together  experts  to  openly  discuss  their 
viewpoints,  a  better  understanding  of  the  subject  results  (Matheson  and 
Howard,  1968:39).  The  face-to-face  method  may  also  bring  out  the 
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critical  issues.  Information  shared  in  this  way  may  spark  new  ideas  or 
a  better  way  of  understanding  old  thoughts. 


Discussion 

These  techniques  are  just  a  few  of  the  methods  available  for 
combining  expert  opinion  and  there  are  plenty  of  variations  on  all  of 
them.  Before  recommending  the  best  method  for  determining  the  branch 
weights,  an  overview  of  the  forms  of  group  opinion,  the  importance  of 
conditional  independence,  the  use  of  sensitivity  analysis,  and  the 
choice  of  experts  will  be  discussed. 

Categories.  (Hogarth,  1986:145)  There  are  three  classes  of  group 
opinion: 


1)  A  group  nay  want  a  consensus  of  opinion  so  as  to  make  a 
decision  within  the  group.  This  could  be  a  husband  and  wife  trying 
to  decide  which  house  to  buy  and  they  need  a  consensus  of  what  they 
both  would  like  to  have  in  their  home. 

2)  Another  form  of  group  opinion  is  the  group  who  needs  a 
consensus  of  opinion  for  presenting  information  to  the  public,  such 
as  in  the  medical  field. 

3)  And  finally  there  is  the  group  opinion  which  will  be  used 
by  a  third  person  for  assessing  a  particular  problem.  This  is  the 
case  for  the  unit  risk  assessment  tree  where  the  third  person  is  a 
decision  maker  who  needs  a  weight  assigned  to  each  branch  of  the 
tree  for  determining  the  weighted  average  of  the  unit  risk. 

This  third  form  of  group  opinion  lends  itself  well  to  evaluation  by  the 

Bayesian  approach  for  combining  expert  opinion.  As  new  evidence  is 

found  for  a  particular  chemical,  the  decision  maker  could  re-evaluate 

the  weights  and  calculate  a  new  unit  risk. 

Conditional  Independence.  If  Independence  between  experts  can  be 

assumed,  then  the  Bayesian  method  simplifies  to  just  an  arithmetic 

average.  By  using  the  face-to-face  interaction  approach  for  combining 


26 


expert  opinion,  the  experts  could  then  be  assumed  to  have  the  saae 
information  available  in  assigning  their  estimated  weights  since  each 
expert  would  be  exposed  to  other  experts'  Information.  Thus, 
conditional  independence  can  be  assumed  and  the  Bayesian  approach  is 
reduced  to  a  simple  expression.  The  Delphi  method  does  allow  for  some 
interaction,  but  the  amount  of  information  transferred  from  one  expert 
to  another  is  limited.  The  face-to-face  allows  for  more  open 
communication  channels  between  experts,  and  a  better  chance  of  each 
expert  understanding  the  others'  reasoning. 

The  differences  of  opinion  left  between  the  experts  after  a  face- 
to-face  debate  could  then  be  averaged  into  a  single  group  opinion  by 
either  arithmetic,  geometric,  or  ranked  (weighted)  averaging  techniques. 
With  the  arithmetic  averaging  method,  if  the  experts  have  a  large  spread 
between  opinion  values,  the  larger  values  will  be  favored.  Since  the 
branch  weights  are  between  zero  and  one,  the  geometric  approach  is  not 
necessary.  Additionally,  "geometric  means  almost  always  emerge  when 
odds  ratios  or  densities  are  used"  and  arithmetic  averaging  is  usually 
used  when  opinions  are  given  in  probability  form  (Cenest  and  Zidek, 
1986:126),  as  is  the  case  here  for  weights.  The  ranked  (weighted) 
averaging  technique  requires  that  the  experts  be  ranked  according  to  the 
decision  maker's  confidence  in  the  experts  ability.  The  decision  maker 
may  do  this  himself  or  may  have  the  experts  rank  themselves.  Assuming 
the  experts  to  be  completely  familiar  with  the  chemical  being  assessed, 
then  ranking  them  may  prove  difficult  and  possibly  embarrassing.  The 
decision  maker  will  most  likely  rank  the  experts  as  equals  which  reduces 
the  ranked  (weighted)  equation  to  the  arithmetic  averaging  technique. 
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Hence,  the  arithmetic  averaging  technique  is  the  regaining  choice  for 

developing  a  combined  group  weight  for  the  tree  branches. 

Sensitivity  Analysis.  Since  it  is  likely  that  the  group  of  experts 

will  not  cone  to  a  complete  agreement  as  to  what  the  weight  for  a 

particular  branch  must  be,  a  decision  maker  should  evaluate  what  effect 

this  disagreement  has  on  the  outcome  of  the  unit  risk.  A  sensitivity 

analysis  could  be  performed  to  assess  the  significance  of  the  experts' 

disagreement  on  the  consensus  (Hogarth,  1975:283). 

Choosing  the  Experts.  Experts  should  be  chosen  who  have  both 

substantive  and  normative  goodness: 

"Substantive  goodness  refers  to  knowledge  which  the  assessor 
[expert]  has  concerning  the  subject  matter  of  concern.  Normative 
goodness  is  the  ability  of  the  assessor  to  express  his  opinions  in 
a  probabilistic  form  [Hogarth,  1975:272]." 

Since  the  choice  of  experts  will  reflect  the  confidence  placed  in  the 
weighted  average  of  unit  risk,  experts  should  be  chosen  on  their  under¬ 
standing  of  the  risk  assessment  procedure  as  well  as  their  knowledge  of 
the  chemical  in  question.  Also,  each  expert  should  have  the  ability  to 
translate  his  preference  for  each  assumption,  to  correctly  estimate  the 
unit  risk  number,  into  a  weight  ranging  from  0  to  1. 

Recommended  Procedure 

From  the  above  discussion,  the  method  for  assessing  the  tree  branch 
weights  is  as  follows: 

1)  Choose  a  panel  of  experts  who  are  familiar  with  the 
chemical  being  assessed  and  have  the  ability  to  translate  their 
beliefs  (preferences)  into  probabilities  (weights). 

2)  Ask  the  experts  to  evaluate  the  tree  diagram  and  assign 
weights  to  the  individual  branches. 
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3)  The  weights  are  then  debated  between  experts  in  a  face-to- 
face  discussion  so  as  to  allow  for  maximum  possible  information 
exchange . 

4)  Any  differences  in  weights  are  then  averaged  using  the 
arithmetic  averaging  technique. 

5)  Finally,  perfora  a  sensitivity  analysis  on  the  tree  branch 
weights,  if  a  consensus  has  not  been  reached,  to  determine  the 
significance  of  the  nonconsensus. 

If  this  procedure  is  followed,  there  will  be  consistency  between 
the  old  assessment  and  a  new  assessment  for  a  particular  chemical  when 
new  information  is  uncovered.  Also,  there  will  be  consistency  from  one 
chemical  to  another.  Consistency  is  required  to  gain  confidence  of  the 
constituencies  who  use  the  unit  risk  number  of  a  chemical  for  evaluating 
hazardous  conditions. 

This  recommended  procedure  for  combining  expert  opinion  can  be  used 
for  combining  the  tree  branch  weights.  In  the  next  chapter  the  dioxin 
and  DCM  tree  diagrams  are  provided  to  a  single  expert  for  assigning  his 
preference  and  likelihood  weights  for  each  assumption  in  the  trees. 

Then  the  trees  are  folded  back  to  piovide  a  weighted  overage  of  unit 
risk  for  DCM  and  a  weighted  average  of  risk-specific  dose  for  dioxin. 
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IV.  PCM  Unit  Risk  Analysis 


In  this  chapter  the  chemical  methylene  chloride  (DCN)  will  have  its 
unit  risk  number  evaluated  using  the  decomposition  method.  Preference 
weights  will  be  assigned  and  a  sensitivity  analysis  performed  to 
determine  the  significant  assumptions. 

Assigning  Preference  Weights 

Dr.  Melvin  Andersen  of  the  Armstrong  Aerospace  Medical  Research 
Laboratory  was  chosen  as  the  expert  for  assigning  preference  weights  to 
each  assumption.  Due  to  time  constraints  and  lack  of  another  convenient 
DCM  expert,  only  one  expert  was  consulted  for  assigning  weights.  Thus 
combining  expert  opinion  as  discussed  in  chapter  3  will  not  be 
illustrated.  The  two  DCM  tree  diagrams  (Figures  3  and  4)  were  provided 
to  Dr.  Andersen  without  unit  risk  numbers.  Dr.  Andersen  then  assigned 
preference  weights  to  each  branch  of  the  tree.  These  weights  are  based 
on  his  belief  of  the  correctness  of  each  assumption  in  estimating  the 
unit  risk  number.  The  preference  weights  for  the  DCM  tree  (Figure  3 
only)  are  listed  in  Table  1. 

According  to  the  preference  weights  assigned  by  Andersen  the 
pharmacodynamic  assumption  weights  are  dependent  on  the  souse 
pharmacokinetic  assumptions.  Also  the  human  pharmacokinetic  assumption 
weights  (Table  1,  column  3)  are  dependent  on  both  the  mouse 
pharmacokinetic  and  species  pharmacodynamic  assumptions,  indicating  that 
depending  on  what  assumption  was  used  for  the  mouse  pharmacokinetics 
does  influence  the  expert's  (Andersen's)  belief  (weight)  of  the  other 
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assumptions  (species  pharmacodynamics  and  human  pharmacokinetics).  The 
preference  weights  for  the  PBPK  (mouse  and  human)  portion  of  the  DCN 
tree  (Figure  4)  are  shown  in  Table  2.  Both  the  mouse  and  human  PBPK 
trees  were  assigned  identical  weights  and  so  only  the  assumptions 
(variables)  with  their  preference  weights  are  shown. 


Table  1.  DCN  Preference  Weights  given  by  Andersen  (Andersen,  1988). 


House 

Species 

Human 

Pharmacokinetics 

Pharmacodynamics 

Pharmacokinetics 

(assumption/ weight) 

(assumption/weight) 

( assumpt i on/we i gh  t ) 

_ > 

PBPK/0.7 

- > 

Body  surface/0.3  — 

Applied/0.3 

->  Applied/0.2  -- 

1 

1 

1 

l 

- > 

Body  weight/0.7  — 

- > 

PBPK/0.7 

! 

i 

» 

- > 

Applied/0.3 

i 

i 

* 

- > 

Body  surface/0.2  --- 

PBPK/1.0 

* 

* 

- > 

Applied/0.0 

->  PBPK/0.8  — 

- > 

PBPK/1.0 

- > 

Body  weight/0.8  — 

- > 

Applied/0.0 

Table  2. 

DCN  Preference  Weights 

for  PBPK  Tree 

(Andersen,  1988). 

Pathway 

Activity  Ratio 

1st  Order 

Maximum  Velocity 

Lung/Liver 

MFO  Rate 

of  Metabolism  by  MFO 

( assumpt 1  on/weight ) 

DCH/0.1 

A1  high/0.5 

Kf  high/0.5 

Vmax  hi*h/0'5 

MFO/0 . 2 

A1  low/0.5 

Kf  low/0.5 

V.ax  low/0’5 

GST/0.7 

A2  high/0.5 

A2  low/0.5 


PBPK  Variable  Analysis 

PBPK  Variable  Values.  The  PBPK  variable  values  for  Al,  A2,  Kf,  and 
V  were  provided  by  Or.  Andersen  <1988),  and  are  shown  in  Table  3. 

ttdX 


Table  3.  PBPK  Variables’  Range  (Andersen,  1988). 


Species 

!  Variable  (units) 

!  Low 

Homlnal 

High 

!  A1 

(ratio) 

!  0.276 

0.405 

0.560 

Mouse 

i  A2 

(ratio) 

!  0.196 

0.280 

0.377 

:  Rf 

(1/hr) 

!  3.000 

4.000 

5.000 

:  vf 

max 

(mg/hr) 

:  0.900 

1.054 

1.200 

!  A1 

(ratio) 

;  o.ooo 

0.001 

0.004 

Human 

:  A2 

(ratio) 

!  0.148 

0.180 

0.297 

:  k. 

(1/hr) 

1  0.320 

0.430 

0.640 

:  vf 

max 

(mg/hr) 

i  85.000 

118.900 

150.000 

PBPK  Sensitivity  Analysis.  The  DCN  tree  as  shown  in  Figures  3  and 
4  create  6546  combinations  for  obtaining  the  unit  risk  number.  Of  the 
6546  combinations  of  assumptions  for  determining  the  unit  risk  number, 
all  but  two  (6544)  are  combinations  using  PBPK  models.  If  any  of  the 
PBPK  variables  could  be  set  to  their  nominal  values  (fixed)  this  would 
reduce  the  number  of  unit  risk  number  calculations  required.  So  a 
stochastic  sensitivity  analysis  (Tatman,  1988)  of  the  PBPK  variables  in 
Table  3  was  performed  in  order  to  determine  the  least  important 
variables  and  hopefully  decrease  the  number  of  variables  in  the  PBPK 
tree.  The  stochastic  sensitivity  analysis  is  described  as  follows: 

The  probabilistic  [stochastic]  sensitivity  indicates  how  the 
certain  equivalent  [weighted  average  of  unit  risk]  depends  on  a 
particular  state  variable  (assumption]  when  the  other  state 
variables  [assumptions]  are  taken  with  their  assigned  probability 
distributions  [weights],  (von  Holstein,  1983:147) 
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The  PBPK  model  developed  by  Andersen  and  Clewell  was  used  to 
generate  the  dose  surrogate  values  for  both  species:  mouse  and  human 
(refer  to  appendix  A  for  details).  The  mouse  dose  surrogate  values  were 
used  in  the  LMS  dose  response  model  (see  appendix  B),  along  with  the  HTP 
bioassay  results  for  mice,  to  obtain  the  dose  response  values.  Then 
these  dose  response  values  were  multiplied  by  the  human  dose  surrogate 
values  (target  dose)  to  get  the  unit  risk  numbers.  Pharmacodynamic 
conversion  (body  surface)  was  not  performed  since  only  the  sensitivity 
of  the  unit  risk  number  to  the  PBPK  tree  variables  was  of  Interest.  Any 
conversion  would  have  been  the  same  for  all  PBPK  variables. 

The  mouse  PBPK  variables  Al,  A2,  K-,  and  V  „  were  combined  with 

the  human  PBPK  variables  Al,  A2,  K, ,  and  V  to  obtain  one  stochastic 

sensitivity  analysis  of  each  PBPK  variable.  The  range  of  unit  risk 
numbers  for  the  stochastic  sensitivity  analysis  is  shown  in  Figure  5  and 
Table  4.  The  stochastic  sensitivity  analysis  shows  that  the  PBPK 
variables  make  very  little  difference  between  their  high  and  low  values 
on  the  unit  risk  number  as  compared  to  the  pathway  assumptions.  For 
pathway  the  range  between  the  high  and  low  unit  risk  numbers  is 
significant.  The  unit  risk  range  for  pathway  is  397  times  larger  then 
the  unit  risk  range  for  variable  Al  (refer  to  Table  4).  The  unit  risk 
number  is  sensitive  to  which  pathway  assumption  is  used  (MFO  or  CST)  and 
is  not  nearly  as  sensitive  to  the  other  four  variables.  Thus  the  PBPK 
variables  Al,  A2,  Kf ,  and  V  were  set  to  nominal  values  as  listed  in 
Table  3. 
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Figure  5  PBPK  Stochastic  Sensitivity  Analysis. 


Table  4. 

Stochastic  Sensitivity  of 

the  PBPK  Variables. 

Variable 

Lower 

Upper  7 

Difference 

Ratio 

unit  risk 

numbers  (xlO  ) 

Pathway 

0.491 

8.420 

7.93 

397 

A2 

0.982 

1.230 

0.248 

12 

Kf 

1.020 

1.180 

0.160 

8 

’max 

1.020 

1.180 

0.160 

8 

A1 

1.090 

1.110 

0,020 

1 

weighted  average  of  unit  risk  *  1.1x10 
for  a  continuous  exposure  to  1  ug/mA3  of  DCM 


Weighted  Average  of  Unit  Risk 

With  tha  four ~PBPK  variables -removed  from  the  tree,  it  reduces  to 
that  shown  in  Figure  6.  How  there  are  only  14  combinations  for 
obtaining  DCH's  unit  risk  number.  In  reality  there  are  28  combinations 
because  a  unit  risk  number  is  calculated  for  both  liver  and  lung  tumors, 
but  these  are  considered  independent  and  are  added  together  to  get  the 
unit  risk  number  for  each  endpoint  of  the  tree. 

The  unit  risk  numbers  associated  with  each  combination  of 
assumptions  in  Figure  6  have  been  determined  (refer  to  appendix  C)  and 
are  listed  in  Table  5.  The  weighted  average  of  unit  risk  is  calculated 
by  converting  the  unit  risk  numbers  for  each  combination  to  logarithms 
and  folding  back  the  tree  using  the  preference  weights  to  obtain  the 
marginal  unit  risk  numbers,  and  the  summation  of  the  marginal  unit  risk 
numbers  is  -6.6748,  which  when  converted  out  of  logarithm  form  gives  an 
unit  risk  number  of  2.1xlO~7. 
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Species 

Pharmaco¬ 

kinetics 


Species  lo 
Hun  jui 
Pharmaco¬ 
dynamics 


Body  Surface 
03  ' 


Applied 

02 


Body  Weight 
07  1 


of  unit  risk  =2.1  E-7 


Body  Surface 
02  ' 


PBPK 

08 


Body  Weight 
08  ' 


Human 

Pharmaco¬ 

kinetics 


PBPK 


Applied 

0.3 


PBPK 

07 

Applied 

0.3 

PBPK 

1.0 

Applied 

0.0 

PBPK 

1.0 


Applied 


Pathway 


0.7 

DCM 


0.7 

DCM 


0.7 

DCM 


0.7 

DCM 


Figure  6.  Reduced  Metbyene  Chloride  Tree  Diagram. 


Unit  Risk 
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1.45  E -6 


1.45  E  -6 


3.48  E -6 


4.46  E -6 


1.15  E-7 


1.15  E  -7 


2.75  E  -7 


1.06  E -5 


555  E  -7 


6.38  E -6 


8.4  E-7 


4.38  E  -8 


5.04  E-7 
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Table  S.  Unit  Risk  Numbers  for  the  DCK  Tree. 


Mouse 

Pharmaco¬ 

kinetics 

Species 
-  Pharmaco¬ 
dynamics 

Human 

Pharmaco¬ 

kinetics 

Pathway 

Unit 

Risk 

Number 

Log 

Unit 

Risk 

Marginal 

Log 

Unit  risk 

applied  body  surface 

PBPK 

MFO 

5.64E-5 

-4.249 

-0.0357 

•« 

i* 

tt 

GST 

1.45E-6 

-5.839 

-0.1717 

•1 

H 

tt 

DCN 

1.45E-6 

-5.839 

-0.0245 

it 

tt 

applied 

n/a 

3.48E-6 

-5.458 

-0.0983 

applied  body  weight 

PBPK 

MFO 

4.46E-6 

-5.351 

-0.1049 

tt 

H 

tt 

GST 

1.15E-7 

-6.939 

-0.4760 

tt 

II 

it 

DCN 

1.15E-7 

-6.939 

-0.0680 

tt 

tt 

applied 

n/a 

2.75E-7 

-6.561 

-0.2755 

PBPK 

body  surface 

PBPK 

MFO 

1.06E-S 

-4.975 

-0.1592 

tt 

tt 

tt 

GST 

5.55E-7 

-6.256 

-0.7006 

It 

it 

tt 

DCM 

6.38E-6 

-5.195 

-0.0831 

It 

tt 

applied 

n/a 

n/a 

n/a 

n/a 

PBPK 

body  weight 

PBPK 

MFO 

8.40E-7 

-6.076 

-0.7777 

tt 

tt 

tt 

GST 

4.38E-8 

-7.359 

-3.2966 

H 

tt 

M 

DCM 

5.04B-7 

-6.290 

-0.4030 

tt 

tt 

applied 

n/a 

n/a 

n/a 

n/a 

log  unit  risk  * 

-6.6748 

UNIT  RISK  • 

2.11E-7 

Using  the  deconposition  method  the  unit  risk  number  for  a 
continuous  exposure  to  1  ug/m3  of  DCN  is  2.1xl0~^.  The  EPA  set  their 
unit  risk  number  to  4.7x10  ^  for  the  same  exposure.  The  decomposition 
method  illustrated  actually  provides  a  slightly  less  conservative  value 
than  does  the  EPA's.  The  DCN  unit  risk  via  the  decomposition  method  is 
approximately  2  times  smaller  than  the  EPA’s  unit  risk. 

Unit  Rl30l.stjrlfr.aUoD 

The  unit  risk  distribution  for  DCN  is  shown  in  Figure  7.  The 
weighted  average  is  dominated  by  the  weights  assigned  to  the  PBPR-body 
weight-GST  assumption  combination.  Almost  45  percent  of  the  weighted 


Unit  Risk  (zlO  ) 


average  Is  Influenced  by  this  combination.  If  the  decision  maker 
believes  the  PBPK-body  weight-CST  combination  should  not  have  this  much 
influence  over  the  weighted  average,  than  he  is  free  to  re-evaluate  the 
unit  risk  number  using  new  preference  weights.  In  this  way  the  decision 
maker  can  easily  analyze  his  unit  risk  number  and  how  it  was  obtained. 

Stochastic  Sensitivity 

Like  the  variables  in  the  PBPK  tree,  the  DCN  tree  can  have  a 
sensitivity  analysis  performed  on  its  assumptions  to  observe  how  much 
they  affect  the  weighted  average.  A  stochastic  sensitivity  analysis  was 
performed  on  the  DCM  tree  and  the  results  are  shown  in  Table  6  and 
Figure  8. 


Table  6.  Stochastic  Sensitivity  of  the  DCM  Tree  Assumptions. 


Assumption 

Lower 
unit  risk 

Upper 

numbers  (xlO  °) 

Difference 

Ratio 

Pathway 

0.0979 

1.73 

1.632 

4.46 

Pharmaco¬ 

dynamics 

0.119 

1.63 

1.511 

4.13 

Human 

Pharmaco¬ 

kinetics 

0.198 

0.589 

0.391 

1.07 

Mouse 

Pharma  CO"- 
kinetics 

0.168 

0.533 

0.366 

1.00 

weighted  average  of  unit  risk  *  2.1E-7 
for  a  continuous  exposure  to  1  ug/m~3  of  DCM 


Of  the  four  assumptions,  pathway  still  influences  the  weighted 
average  more  than  the  other  three,  but  not  as  significantly  as  it  did  in 
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20.0- 


e  8.  DCM  Unit  Risk  Sensitivity  Analysis. 


the  PBPK  tree  where  it  significantly  dominated  the  four  PBPR  variables 
(4.46  tines  larger  range  as  compared  to  397  times  larger  range).  Hence, 
the  four  assumptions,  pathway,  pharmacodynamics,  human  pharmacokinetics, 
and  mouse  pharmacokinetics,  are  important  in  the  unit  risk  assessment 
process,  at  least  with  DCM. 

Discussion 

The  unit  risk  number  of  DCM  for  a  continuous  exposure  to  1  ug/m^ 
was  found  to  be  2.1xl0~7  via  the  decomposition  method.  Even  though  the 
2.1xl0~7  seems  like  a  good  unit  risk  number  it  is  not  a  recommendation 
of  the  Air  Force,  but  rather  an  example  to  illustrate  the  decomposition 
method  for  determining  the  risk  associated  with  chemicals.  The 
significant  assumptions  for  determining  the  unit  risk  in  order  of 
decreasing  significance  are:  the  PBPK  pathway,  pharmacodynamics,  human 
pharmacokinetics,  and  mouse  pharmacokinetics.  The  PBPK  variables  Al, 

A2,  K.,  and  V  do  not  make  much  difference.  This  study  indicates  that 

t  max 

scientific  efforts  should  concentrate  their  efforts  on  analyzing  the 
significant  assumptions  since  the  uncertainty  of  these  play  a  larger 
role  in  the  variance  of  the  unit  risk  number. 

The  decomposition  method  can  be  applied  when  evaluating  the  unit 
risk  number  of  DCM.  This  method  incorporates  all  acceptable  methods  for 
determining  the  unit  risk  of  DCM  whereas  the  EPA's  unit  risk  is 
developed  from  only  one  combination  of  these  assumptions.  The 
decomposition  method  ranks  each  assumption  as  to  its  ability  to  predict 
the  unit  risk  of  the  chemical  while  incorporating  the  most  current 
information  available  to  the  toxicologist. 
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Dioxin  is  used  to  illustrate  the  flexibility  of  the  decomposition 
method.  Instead  of  calculating  a  weighted  average  of  unit  risk,  a 
weighted  average  of  risk-specific  dose  (RSD)  will  be  determined  for 
dioxin.  Unlike  the  unit  risk,  which  is  a  risk  associated  with  a  set 
exposure,  the  risk-specific  dose  is  an  exposure  associated  with  a  set 
risk.  The  dioxin  tree  (Figure  2)  will  have  weights  assigned  to  its 
branches  and  a  weighted  average  of  risk-specific  dose  calculated  for  a 
set  risk  probability  of  lxl0~6. 


Assigning  Preference  Weights 

Again  Dr.  Andersen  was  chosen  as  the  expert  for  assigning 
preference  weights  to  each  assumption.  The  dioxin  tree  diagram, 

Figure  2,  was  provided  to  Dr.  Andersen;  he  then  assigned  his  preference 
weights  to  each  branch,  which  are  shown  in  Figure  9. 


Weighted  Average  of  Risk-Specific  Dose 

There  are  14  combinations  for  obtaining  weighted  average  of  risk- 
specific  dose  for  dioxin  (EPA,  November  1987).  The  weighted  average  of 
risk-specific  dose  can  now  be  calculated  using  the  RSD  associated  with 
each  combination  of  assumptions  with  their  corresponding  preference 
weight  in  Figure  9.  The  weighted  average  is  calculated  by  converting 
the  RSDs  for  each  combination  to  logarithms  and  folding  back  the  tree 
using  the  preference  weights. 

Using  Dr.  Andersen's  preference  weights,  the  decomposition  method 
places  the  risk-specific  dose  for  dioxin  at  0.7  pg/kg/da y.  Which  means 
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Risk-specific 


0.7  pg/kg/day 


Risk-Specific 
Dose  (pg/kg/day) 

Safety  Factor 

100  HAW 

ra -  100 
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05 
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05 
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Figure?,  Dioxin  Tree  Diagram  w/ Preference  Weights. 
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that  for  a  health  risk  of  1x10  6  a  person  nay  receive  an  internal 
exposure  of  0.7  pg/kg/day.  This  converts  to  an  external  exposure  doss 
of  2.45xl0~6  ug/*3  using  the  applied  dose  conversion.  The  EPA  set  their 
risk-specific  dose  at  0.1  pg/kg/day  for  the  same  health  risk. 

Discussion 

The  decomposition  method  illustrated  actually  provides  a  slightly 
less  conservative  value  than  does  the  EPA's.  The  dioxin  RSD  via  the 
decomposition  method  is  approximately  7  times  higher  than  the  EPA's 
allowed  dose.  Of  course  different  preference  weights  would  generate  a 
different  risk-specific  dose. 

Again  this  calculated  value  for  risk-specific  dose  is  not  a  recom¬ 
mendation  of  the  Air  Force,  but  just  an  example  to  illustrate  the  decom¬ 
position  method  for  determining  the  risk  associated  with  chemicals.  The 
decomposition  method  can  be  applied  for  evaluating  the  risk-specific 
dose  of  dioxin.  This  method  incorporates  all  acceptable  methods  for 
determining  the  risk-specific  dose  of  dioxin  whereas  the  EPA's  risk- 
specific  dose  was  developed  from  a  subjective  debate  of  the  other 
methods  and  a  conservative  middle  ground  number  of  0.1  pg/kg/day  was 
chosen.  The  decomposition  method  ranks  each  assumption  as  to  its 
ability  to  predict  the  risk  specific  dose  of  the  chemical  while  incorpo¬ 
rating  the  most  current  information  available  to  the  toxicologist. 

When  pharmacokinetic  models  are  developed  for  dioxin,  the  dioxin 
tree  could  easily  be  modified  to  include  them.  The  flexibility  of  the 
decomposition  method  lends  itself  for  many  different  risk  analysis 
methods  and  allows  updating  as  new  information  unfolds. 
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VI.  Suraary.  Conclusions  and  Recommendations 


Summary 

There  are  various  ways  of  deteraining  the  risk  of  a  cheaical.  The 
different  federal  agencies  responsible  for  setting  cheaical  risk 
standards  are  not  in  coaplete  agreement.  Overestimating  the  risk  places 
a  financial  burden  on  industry  for  maintaining  the  standard,  and 
underestimating  the  risk  puts  society  in  jeopardy  of  a  greater  cancer 
threat.  This  study  has  shown  how  the  different  risk  assessment  methods 
can  be  decomposed  into  their  assumptions;  these  assumptions  are  then 
evaluated  for  their  ability  to  accurately  estimate  the  risk  of  a 
chemical.  Then  the  weighted  assumptions  are  recomposed  into  a  weighted 
average  of  risk. 

Conclusions 

The  decomposition  method  can  be  applied  to  evaluating  the  unit  risk 
of  chemicals.  This  method  incorporates  several  acceptable  unit  risk 
methods  into  one  larger  and  more  thorough  process  instead  of  just 
relying  on  one  method,  which  may  ignore  the  uncertainty  associated  with 
its  assumptions.  In  doing  so,  the  significance  of  the  assumptions  and 
variables  which  are  used  can  be  appraised  and  compared.  Each  assumption 
is  weighted  according  to  its  ability  of  estimating  the  unit  risk  number 
and  through  this  weighting  an  objective  assessment  of  the  chemical's 
risk  is  determined.  Finally,  a  weighted  distribution  of  the  unit  risk 
number  can  be  generated  for  providing  a  better  understanding  of  how  the 
unit  risk  number  was  obtained  and  the  possible  range  of  values. 
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The  decomposition  method  aids  the  decision  maker  in  better 
understanding  the  chemical  in  question.  By  addressing  the  various 
assumptions  used  for  determining  the  unit  risk  number,  the  decomposition 
method  forces  the  decision  maker  to  place  a  likely  range  of  values  and 
preference  weights  on  the  assumptions.  In  doing  this  the  decision  maker 
deals  with  the  uncertainty  of  each  assumption  with  the  hope  of  obtaining 
a  realistic  and  acceptable  unit  risk  number. 

Recommendations 

When  toxicologist  develop  a  physiologically-based  pharmacodynamic 
model,  replacing  body  weight  and  body  surface  methods,  then  the 
decomposition  method  should  be  used  to  evaluate  the  significance  of  the 
new  model  and  generate  a  new  unit  risk  for  DCM,  or  any  other  chemical 
that  has  had  a  PBPK  model  built  for  it. 

The  pathway  the  chemical  < 1 . e . ,  DCM)  uses  to  attack  the  body  should 
be  investigated  further.  By  reducing  the  uncertainty  associated  with 
the  pathway,  the  unit  risk  distribution  of  a  chemical  could  be  greatly 
narrowed.  The  same  goes  for  pharmacokinetics  in  human  and  laboratory 
animals. 

The  mixing  of  preference  and  likelihood  weights  (probabilities) 
needs  to  be  addressed.  The  recomposition  of  the  weighted  assumptions 
and  probability  distributions  of  the  PBPK  variables  into  one  chemical 
ri3k  poses  theoretical  obstacles.  Future  research  into  this  area  of 
mixing  weights  and  probabilities  is  required  so  as  to  provide  a 
theoretically  sound  unit  risk  number.  As  for  now,  the  decomposition 
method  splits  the  tree  diagram  into  two  areas:  preference  weights  and 
probabilities. 
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Appendix  A:  The  Physiologically-Based  Pharmacokinetic  Model 

for  Methylene  Chloride 


The  physiologically-based  pharaacoklnetlc  (PBPK)  model  for 
■ethelyene  chloride  (DCN)  was  developed  by  Clewell  and  Andersen  of  the 
Armstrong  Aerospace  Medical  Research  Laboratory  at  Wright  Patterson  AFB, 
Ohio  (Andersen  et  al.,  1987).  The  parameter  values  used  for  both  mouse 
and  human  are  as  follows: 


MOUSECS 

"MOUSE  DOSE  SURROGATES") 

MOUSE) 

SET  BW  =  0 .0345) 

SET  QCC*14.4,  QPC=14.4) 

SET  VLC=0 . 04 , VRC=0 . 05 ,VFC=Q .04 ,VSC=0. 78) 

SET  QLC=0 . 24 ,QFC=0 . 05 ,QRC=0 . 52 ,QSC=0 .19) 

SET  PL=1 . 71 , PF=14.5 , PS=0 .960 , PR=1 . 71 , PLU=1 . 71 , PB=8 . 29) 

SET  AU0.405,  A2=0.28) 

SET  RF=4.0172,  VMAX*1.0543,  KM*0.396) 

SET  TITLE  =  "PB-PK  MODEL  OF  MECL2") 

SET  TCHNG=6 , TSTOP=24 , I VDOSE=0) 

START) 

HUMANDS 

"HUMAN  DOSE  SURROGATES") 

HUMAN) 

SET  BW=70 . 0) 

SET  QCC=24,  QPC  =  24) 

SET  VLC=0.0314,  VFC=0.231,  VRC=0.0371,  VSC=0.621) 

SET  QLC*0.24,  QFC=0.05,  QRC*0.52,  QSC=0.19) 

SET  PB=9 . 7 ,  PL=1 .46,  PR=1.46,  PS=0.82,  PF=12.4,  PLU=1.46) 
SET  Al-0.001,  A2=0 . 18) 

SET  VMAX=118.9,  KM=0.58,  KF-0.43) 

SET  TITLE  =  "PB-PK  MODEL  OF  MECL2") 

SET  TCHNG“6  ,TST'>P=24 ,  IVDOSE*0) 

START) 


Computer  output  for  a  "MOUSEDS"  is  as  follows: 

D  CONC, DRINK) 

CONC  4000.00000  DRINK  0. 

D  RISK1L,RISK1P) 

RISK1L  5474.81000  RISK1P  2369.07000 


D  RISK2L,RISK2P! 

RISK2L  2292.31000  RISK2P  683.985000 

0  AUCL, AUCLU! 

AUCL  975.769000  AUCLU  1039.83000 

D  COHC, DRINKS 

CONC  2000.00000  DRINK  0. 

D  RISK1L.RISK1P! 

RISK1L  5154.93000  RISK1P  2237.63000 

n  Brsr?i 

RISK2L  1001.35000  RISK2P  313.775000 

D  AUCL, AUCLU! 

AUCL  426.243000  AUCLU  477.017000 


Computer  output  for  a  "HUMANDS"  is  as  follows: 

D  CONC, DRINK! 

CONC  4000.00000  DRINK  0. 

D  RISK1L, RISK1P! 

RISK1L  1195.38000  RISK1P  1.21114000 

D  RISK2L,RISK2P! 

RISK2L  297.479000  RISK2P  55.4503000 

D  AUCL, AUCLU! 

AUCL  1010.05000  AUCLU  1045.96000 

D  CONC, DRINK! 

CONC  2.8764E-04  DRINK  0. 

D  RISK1L.RISK1P! 

RISK1L  8 . 3471E-04  RISK1P  2.1215E-06 

D  RISK2L,RISK2P! 

RISK2L  3.8484E-06  RI5K2P  1.7606E-06 

D  AUCL, AUCLU! 

AUCL  1 . 3067E-05  AUCLU  3.3210E-05 


The  computer  command  file  for  running  the  DCM  PBPK  model  (SAB.CSL) 
is  as  follows  (Andersen  et  al.,  1987): 


"SAB.CMD" 

•'Command  File  for  SAB.CSL  -  Rev.  11/18/85” 

"  To  Generate  Data  for  SAB  Risk  Analysis  " 

"  Can  be  used  for  inhal  if  CONC  >  0,  or” 

■'  for  Drinking  water  by  setting  DRINK  to  PPM  in  H20" 

PREPAR  T , CA , CV , CA1 , CX , RALU , ALU , CLU , AUCLU 
PREPAR  RAMI P , AMI P , RAM2P , AM2 P , CF , CL , AUCL , CS , CR 
PREPAR  RAM 1 L , AMI L , RAM2L , AM2L , CV , CVF , CVL , CVS , CVR 
SET  WESITC=.F. ,FTSPLT=.T. 
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PROCED  HOUSE 

SET  BW  =  0.0345 

SET  QCC=14.4,  QPC*14.4 

SET  VLC-0 . 04 , VRC=0 . 05 , VFC=0 . 04 , VSC=Q . 78 

SET  QLC=0 . 24 ,QFC=0 . 05 ,QRC=0 . 5  2 ,QSC=0 . 19 

SET  PL=1.71,PF=14.5,PS=0.960,PR=1.71,PLU=1.71,PB=8.29 

SET  A1=0 .405 ,  A2=0.28 

SET  KF=4.0172,  VMAX=1.0543,  KM=0.396 

END  $"END  OF  PROCED  B6C3F1-K0USEH 

PROCED  HUNAN 

SET  BW*70.0 

SET  QCC=24,  QPC  =  24 

SET  VLC=0.0314,  VFC=0.231,  VRC=0.0371,  VSC=0.621 

SET  QLC=0.24,  QFC=0.05,  QRC=0.52,  QSC*0.19 

SET  PB=9.7,  PL=1 .46 ,  PR=1.46,  PS=0.82,  PF=12.4,  PLJ=i.46 

SET  Al=0.001,  A2=0 . 18 

SET  VMAX-118.9,  KM=0.58,  KF=0.43 

END  I'END  OF  PROC  HUMAN' 

PROCED  HOUSEV  I* MOUSE  VALIDATION' 

MOUSE 

SET  TITLE  =' MOUSE  —  50  &  10  MG/KG  IV 
SET  CONC=0 , I VDOSE=50 , TSTOP* 1.5, TCHNG=6 . 0 
START 

SET  NRWITGs.T. , IVDOSE=10 
START 

PLOT  CV, 'LOG* , 'LO' =.001 , 'HI '=100, 'XHI '=1.5 
SET  NRUITGs . F . , I VDOSE=0 
END  $ 'END  OF  PROC  MOUSEV' 

PROCED  HUMANV  $' HUMAN  VALIDATION’ 

HUMAN 

SET  TITLE  =  'HUMAN  —  350  &  100  PPM  INH 

SET  CONC=350,TCHNG*6 ,TSTOP=36 

START 

SET  CONC=100,NRWITG=.T. 

START 

PLOT  CV, 'LOG' , 'XHI ' =36 , ' LO' = . 01 , 'HI ' =100 
SET  NRWITG-.F. 

END  I'END  OF  PROC  HUMANV' 

PROCED  MOUSEDS  | MOUSE  DOSE  SURROGATES' 

MOUSE 

SET  TITLE  *  ' PB-PK  MODEL  OF  MECL2' 

SET  TCHNG=6 ,TSTOP=24, IVDOSE=0 
START 

D  CONC, DRINK 
D  RISK1L.RISK1P 
D  RISK2L,RISK2P 
D  AUCL.AUCLU 

END  I'END  OF  PROC  MOUSEDS' 
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PROCED  HUMANDS  $' HUMAN  DOSE  SURROGATES* 

HUNAN 

SET  TITLE  =  ' PB-PK  MODEL  OF  MECL2  * 

SET  TCHNC»6,TSTOP«24,IVDOSE«0 
START 

D  CONC, DRINK 
D  RISK1L.RISK1P 
D  RISK2L,RISK2P 
D  AUCL.AUCLU 

END  I 'END  OF  PROC  HUMANDS* 

PROCED  DS 
D  CONC, DR INK 
D  RISK1L.RISK1P 
D  RISK2L, RISK2P 
D  AUCL,AUCLU 
END  I'END  OF  PROC  DS* 

PROCED  MENU 

t  • 

•ENTER  ONE  OF  THE  FOLLOWING:* 

I  • 

'MOUSEV  —  plot  mouse  validation  runs' 

'HUMANV  —  plot  human  validation  runs' 

•  • 

'MOUSEDS  —  calculate  mouse  dose  surrogates  (SET  CONC  and 
DRINK  first)* 

'HUMANDS  —  calculate  human  dose  surrogates  (SET  CONC  and 
DRINK  first)’ 

' DS  —  display  dose  surrogates' 

'MENU  --  type  this  menu' 

•  • 

END  I'END  OF  PROC  HELP’ 

SET  HVDPRN* . T . , GRDCPL= . F . , C JV ITG= . F . 

MENU 

SET  CMD«*5 
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The  PBPK  computer  model  (SAB.CSL)  is  as  follows  (Andersen  et  al.,  1987) 


PROGRAM  SAB.CSL  -  Revision  =  11/18/05  at  02:21  pm 

"This  program  will  handle  inhalation  or  drinking  water  Exposures 
"To  use  inhalation,  set  the  variable  COHC  to  PPM,  DRINK  to  0.0" 
"To  use  drinking  water,  set  DRINK  to  PPM  in  water,  CONC  to  0.0" 
"Default  is  inhalation  exposure/B6C3Fl  mouse/4000  PPM" 

"set  CHANGE  =  .F.  to  maintain  cint  constant  during  runs" 

"Will  handle  iv  doses  if  you  set  IVDOSE  >  0.0" 

INITIAL 

LOGICAL  CHANGE  $"T0  CONTROL  ALTERATION  OF  CINT  DURING 
RUNS" 

CONSTANT  TINF  =0.01 
CONSTANT  IVDOSE  =0.0 


CONSTANT 

QPC=28 .0 

CONSTANT 

QCC=20.O 

CONSTANT 

VLC=0.04 

CONSTANT 

VFC=0.04 

CONSTANT 

VRC=0.05 

CONSTANT 

VSC=0. 78 

CONSTANT 

VLUC=0 .0115 

CONSTANT 

QLC=0 .24 

CONSTANT 

QFC=0 .05 

CONSTANT 

QSC=0. 19 

CONSTANT 

QRC=0 .52 

CONSTANT 

PL  =  1.71 

CONSTANT 

PF  =  14.5 

CONSTANT 

PS  *  0.96 

CONSTANT 

PLU  =1.71 

CONSTANT 

PR  =  1.71 

CONSTANT 

PB  =  8.29 

"Constants  for  Mouse" 

("Unsealed  Alveolar  Vent" 
("Unsealed  Cardiac  Output" 

("Vol  Liver  as  X  Body  Wt" 

("Vol  Fat  as  X  Body  Wt" 

("Vol  Rapid  Per.  as  X  Body  Wt" 
("Vol  Slow  Per.  as  X  Body  Wt" 
("Unsealed  vol  of  lung  tissue" 

("Flow  to  Liver  as  X  Cardiac  Out.” 
("Flow  to  Fat  as  X  Cardiac  Out." 
("Flow  to  Slow  as  X  Cardiac  Out." 
("Flow  to  Rapid  as  X  Cardiac  Out." 

("Liver/Blood  Part  Coeff." 
("Fat/Blood  Part  Coeff." 
("Slow/Blood  Part  Coeff." 
("Lung/Arterial  Blood  Part  Coeff.” 
("Rapid/Blood  Part  Coeff." 
("Blood/Air  Part  Coeff." 


CONSTANT  CONC  =  4000.0 
CONSTANT  BW  =  0.0345 
CONSTANT  VNAX  =  1.054 
CONSTANT  KM  =  0.396 
CONSTANT  KF  =  4.017 
CONSTANT  Al  =  0.416 
CONSTANT  A2  =  0.137 
CONSTANT  DRINK  =  0.0 
CONSTANT  CHANGE  =  .TRUE. 
CONSTANT  CONC 2  =  0. 


("Cone,  of  MeCL2  in  PPM" 

("Body  Wt  in  kg" 

("Unsealed  VMax  for  Sat  Pathway" 
("Km  for  Sat  Pathway" 

("1st  Order  Rate  Constant" 

("VMax (Lung)  /  VMax(Liver)" 
("KF(Lung)  /  KF( Liver)" 
("Concentration  in  PPM  in  H20" 
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"**  TIMING  COMMANDS" 


CONSTANT  TSTOP  =24.0 

CONSTANT  POINTS  =  96.0 

CONSTANT  TCHNG  =  6.0 

CINT  =  TSTOP  /  POINTS 

"SCALED  PARAMETERS" 

QC  =  QCC*BW**0. 74 
QP  =  QPC*BW**0.74 
QL  =  QLC*QC 
QF  ■  QFC*QC 
QS  =  QSC*QC 
QR  =  QRC*QC 
VL  =  VLC*BW 
VF  =  VFC*BW 
VS  =  VSC*BW 
VR  =  VRC*BW 
VLU  =  VLUC  *  BW**.99 

lUKR  =  uRiNR/24  *  0.102  *  BW**0.7 

IVR  =  IVDOSE  *  BW  /  TINF 


END  $"£ND  OF  INITIAL" 


DYNAMIC 


ALGORITHM  IALG  =  2 
T2  =  T  *  60. 


IF  (.NOT.  CHANGE)  GOTO  NOCH 
CONC3  =  CONC  *  (1.0  -  STEP(TCHNG) ) 

IF  (CONC2  .NE.  CONC3)  CINT  =  0.0001 
IF  (CONC2  .EQ.  CONC 3)  CINT  =  TSTOP  /  POINTS 
CONC2  =  CONC 3 
NOCH..  CONTINUE 

DERIVATIVE 

"Calculate  Concentration  Inhaled" 

Cl  =  0.00348  *  CONC  *  (1.0  -  STEP(TCHNG) )  $"Conc  in  mg/1" 

"Algebraic  Solution  for  CA1  after  gas  exchange" 

CA1  =  (QC*CV  +  QP*CI )  /  (QC  +  QP/PB) 

CX  =  CA1/PB 
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"Mass  Balance  for  the  Lung  Tissue  Compartment" 
RALU  *  QC*(CA1  -CA)  -  RAM1P  -  RAM2P 
ALU  =  IHTEC< RALU, 0.0) 

CLU  =  ALU  /  VLU 
CA  -  CLU  /  PLU 
AUCLU  *  INTEG(CLU,  0.0) 

RAM1P  =  A1  *  VMAX  *  CA  *  VLU  /  VL  /  (KM  ♦  CA) 

AM1P  =  INTEG(RAM1P,0.0) 

RAM2P  =  A2  *  KF  *  CA  *  VLU 
AM2P  *  INTEG(RAM2P,0.0) 

"AF  3  Amount  in  Fat  Compartment  (mg)" 

RAF  =  QFMCA-CVF) 

AF  *  INTEG(RAF,  0.) 

CVF  *  AF/(VF*PF) 

CF  =  AF/VF 

"AL  =  AMOUNT  IN  LIVER  COMPARTMENT  (MG)" 

RAL  =  QL*(CA-CVL)  -  RAMIL  -  RAM2L  *  KZER 
AL  =  INTEG(RAL,  0.) 

CVL  =  AL/( VL*PL) 

CL  =  AL/VL 
AUCL  =  INTEG(CL.O.O) 

“AS  *  AMOUNT  IN  SLOWLY  PERFUSED  TISSUES  (MG)" 
RAS  *  QS*(CA  -  CVS) 

AS  =  INTEG(RAS,  0.) 

CVS  =  AS/(VS*PS) 

CS  a  A S/VS 

"AR  =  AMOUNT  IN  RAPIDLY  PERFUSED  TISSUES  (MG)" 
RAR  =  QR*(CA  -CVR) 

AR  =  INTEG(RAR,  0.) 

CVR  =  AR/( VR*PR) 

CR  =  AR/VR 

"ANIL  &  AM2L  -  Amounts  metabolized  in  LIVER" 
RAMIL  -  VMAX  *  CVL  /  (KM  +  CVL) 

AM1L  -  INTEG( RAMIL, 0.0) 

RAM2L  =»  KF*CVL*VL 
AM2L  *  INTEG(RAM2L,0.0) 

IVZONE  «  RSW(T.GE.TINF,  0.0,  1.0) 

IV  *  IVZONE  *  IVR 
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"CV  -  MIXED  VENOUS  BLOOD  CONCENTRATION  (MG/L)" 
CV  =  (QF*CVF  f  QL*CVL  +  QS*CVS  +  QR*CVR  t  IV)/QC 

TERMT(T.GE.TSTOP) 

END  $"END  OF  DERIVATIVE" 

END  $"END  OF  DYNAMIC" 

TERMINAL 

"RISK  FACTORS  AS  FUNCTION  OF  VOLUME  OF  TISSUE" 

RISK1L  =  AM1L/VL 
RISK2L  =  AM2L/VL 
RISK1P  *  AM1P/VLU 
RISK2P  =  AM2P/VLU 

END  $"END  OF  TERMINAL" 

END  $"END  OF  PROGRAM" 
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Appendix  B:  Linear  Multi-Stage  Dose  Response  Model 


The  linear  multi-stage  (LMS)  model  is  used  by  most  federal  agencies 
for  calculating  the  dose  response.  The  LMS  model  used  for  calculating 
the  dose  response  for  DCM  is  the  GLOBAL83.FOR  program  provided  by  AAMRL 
and  is  partially  discussed  as  follows  (Howe  and  Crump,  1983): 


To  run  GLOBAL83.FOR,  edit  the  file  GLOBAL83.DAT.  The  file  looks 
like  this: 


Methylene  Chloride  Risk  Assessment 

3, -1,0, 1,0, 0,0, 2 

50,48,48 

0,231,482 

3,30,41 


and  the  meaning  of  the  lines  is  given  below: 

C ***************************************************************** 

c 

c  THE  INPUT  IS  AS  FOLLOWS: 

CARD  1: TITLE 

CARD  2 : LEVELS , NDGRE , NOZERO , I ADD , NAX , NAY , I BAY , LU 

CARD  3:NUM( I ) , 1=1 , LEVELS 

CARD  4:DOSE( I ) , 1=1 .LEVELS 

CARD  5 : 1  DEAD(  I ) ,  I  *  1 ,  LEVELS 

C  IF  NAX  NOT  ZERO 

CARD  6 : (AX( I ) , IAX( I ) ) , 1=1 , NAX 

C  IF  NAY  NOT  ZERO 

CARD  7 : (AY( I) , IAX( I ) ) , 1=1 , NAY 

C 

C  ALL  CARDS  ARE  FREE  FORMAT  EXCEPT  TITLE 

C 

C  LEVELS  IS  THE  NUMBER  OF  DOSE  VALUES 
C  NDGRE  IS  THE  DEGREE  OF  POLYNOMIAL  DESIRED 
C  IF  NDGRE  IS  NEGATIVE  (MULTISTAGE  OPTION)  THE  PROGRAM 

C  AUTOMATICALLY  SETS  THE  DEGREE  TO  LEVELS  MINUS  ONE 

C  NOZERO  SHOULD  BE  ONE  IF  IT  IS  KNOWN  THAT  THERE  IS  NO  BACKGROUND 
C  OTHERWISE  IT  SHOULD  BE  ZERO 
C  IADD=1  FOR  ADDITIONAL  RISK**  P(D)-P(0) 

C  IADD=2  FOR  EXTRA  RISK**  ( P( D)-P(0) )/( 1-P(0) ) 

C  NAX  IS  THE  NUMBER  OF  RISKS  FOR  WHICH  MLE  DOSE  AND  DOSE  LOWER 
BOUND 

C  AND/OR  UPPER  BOUND  ARE  DESIRED. 

C  NAY  IS  THE  NUMBER  OF  DOSES  FOR  WHICH  MLE  RISK  AND  UPPER  BOUND  ON 
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C  RISK  IS  DESIRED. 

C  IBAY=0  IF  GLOBAL  79  CONFIDENCE  INTERVALS  ARE  ALSO  DESIRED. 

C  I BAY* 1  IF  THIS  OUTPUT  IS  NOT  WANTED 

C  LU  =0  IF  LOWER  BOUNDS  ON  DOSE  ARE  DESIRED 

C  LU  si  IF  UPPER  BOUNDS  ON  DOSE  ARE  DESIRED 

C  LU  =2  FOR  BOTH  UPPER  AND  LOWER  BOUNDS 

C  NUH(I)  IS  THE  NUMBER  OF  ANIMALS  ON  TEST  IN  I-TH  CROUP 
C  DOSE(I)  IS  THE  DOSE  VALUE  FOR  THE  I-TH  GROUP 
C  IDEAD(I)  IS  THE  NUMBER  OF  DEFECTIVE  ANIMALS  AT  I-TH  LEVEL 
C  AX ( I )  ARE  RISK  AT  WHICH  SAFE  DOSE  LOWER  BOUNDS  ARE  CALCULATED 
C  AY(I)  ARE  DOSES  AT  WHICH  UPPER  BOUNDS  ON  RISK  ARE  CALCULATED 
C  I AX (I)  INDICATES  THE  SIZE  OF  CONFIDENCE  INTERVALS. 

C  IAXsl  FOR  90%  ;IAX=2  FOR  95%  ;IAX=3  FOR  97. 5%  ;IAX=4  FOR  99X 

C 

C***************************************************************** 


When  the  CLOBAL83.DAT  file  is  ready  Just  type  @GL  to  run  the 
GLOBAL83.FOR  program.  The  output  will  be  typed  to  the  screen  and 
stored  in  the  file  GL0BAL83.0UT.  The  unit  risk  can  be  calculated 
froa  the  last  entry  in  the  output,  which  looks  like  this: 


CONFIDENCE  LIMITS  FOR  A  RISK  OF  0.1E-07  M.L.E.  DOSE*0.2717277724E-05 

UPPER  CONFIDENCE  LIMITS  ON  RISK 

90%  95%  97. 5%  99% 

0 . 120090E-07  0. 126223E-07  0.131714E-07  0.138304E-07 

LOWER  CONFIDENCE  LIMITS  ON  SAFE  DOSE 
90%  95%  97.5%  99% 

0 . 226271E-05  0.215275E-05  0.206302E-05  0.196471E-05 


To  obtain  the  unit  risk,  divide  the  risk,  0.1E-7,  by  the  95%  lower 
confidence  limit  on  the  dose,  0.215E-5,  to  get  4.65E-3.  If  the  dose 
used  in  the  GLOBAL83.DAT  file  is  a  dose  surrogate,  as  in  this  case,  then 
run  the  PBPK  model  for  humans  at  1  ug/m3  and  multiply  the  resulting 
prediction  for  the  dose  surrogate,  6.16E-6,  by  the  risk,  0.1E-7,  and 
divide  by  the  95%  lower  confidence  limit  on  dose,  0.215E-5,  to  get  a 
unit  risk  of  2.86E-8. 
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Appendix  C:  Calculating  PCM  Unit  Risk  Humbers 

The  14  combinations  for  determining  the  unit  risk  number  of  DCM  are 

shown  in  Figure  6.  The  various  conversion  factors  used  in  determining 

the  unit  risk  are  listed  below: 

Human  weight  -  70  kg 

Human  breathing  rate  =  20  m*3/day 

Human  exposure  to  1  ppm  a  0.993314285  mg/kg/day 

House  weight  a  0.0345  kg 

Mouse  rate  =  0.043  sT3/day 

House  exposure  to  1  ppm  *  0.728261113  mg/kg/day 

Body  Surface  conversion  factor  ■  12.65978455 

The  unit  risk  numbers  for  the  14  combinations  are  developed  for 
both  lung  and  liver  tumor  cases  from  the  NTP  bloassay  study  (1986).  The 
unit  risk  numbers  calculated  for  lung  and  liver  are  considered 
independent  and  are  added  together  to  create  one  unit  risk  number  for 
each  of  the  14  combinations  of  assumtions  from  Figure  6.  The  methods 
Implemented  for  calculating  the  unit  risk  for  each  combination  of 
assumptions  were  outlined  by  Dr.  Andersen  and  Lt.  Col.  Harvey  Clewell  of 
AAMRL  and  Dr.  Murray  Cohn  of  the  Health  Sciences  Directorate  in 
Washington  DC.  The  assumption  combinations  are  listed  in  the  following 
order:  mouse  pharmacokinetic  -  species  to  human  pharmacodynamic  -  human 
pharmacokinetic.  The  unit  risk  calculations  are  shown  on  the  following 
pages. 
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Applied  dose  -  Body  weight  (and  Body  surface)  -  Applied  dose 


combinations  (Andersen  and  Clewell,  1988): 


EFFECT 

NTP 

BIOASSAY 

RESPONSE 

MOUSE 

EXPOSURE 

DOSE 

MOUSE 

INTERNAL 

GLOBAL83 

HUMAN 

INTERNAL 

GLOBAL83 

95X  LCL  at 
1.00E-08 

tumor/group 

ppm 

mg/kg/day 

mg/kg/day 

risk 

LUNG 

3/50 

0 

0 

0 

30/48 

2000 

1457 

1457 

1 . 3272E-05 

41/48 

4000 

2913 

2913 

LIVER 

3/50 

0 

0 

0 

16/48 

2000 

1457 

1457 

4 . 8101E-05 

40/48 

4000 

2913 

2913 

EFFECT 

HUMAN 

EXPOSURE 

DOSE 

HUMAN 

INTERNAL 

DOSE 

TARGET 

DOSE 

UNIT  RISK 

95 X  LCL/ 
Target  dose 

COMBINED 

UNIT  RISK 
lung+ liver 

LUNG 

ppm 

2.88E-04 

mg/kg/day 

2.8571E-04 

mg/kg/day 

2.8571E-04 

risk 

2.15E-07 

BODY  WEIGHT 
2.75E-07 

LIVER 

2.88E-04 

2.8571E-04 

2.8571E-04 

5.94E-08 

BODY  SURFACE 
3.48E-06 
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Applied  dose  -  Body  weight  (and  Body  surface)  -  PBPR  method 


combinations  (Cohn,  1987): 


PATHWAY/ 

EFFECT 

NTP 

BIOASSAY 

RESPONSE 

MOUSE 

EXPOSURE 

DOSE 

MOUSE 

DOSE 

SURROGATE 

DOSE  SUR. 
MOD 
FACTOR 

MOUSE 

INTERNAL 

GST/LUNG 

30/48 

ppB 

2000 

mg/ltr/day 

314 

1.09 

mg/kg/day 

1457 

41/48 

4000 

684 

1.00 

2913 

CST/LIVER 

16/48 

2000 

1001 

1.14 

1457 

40/48 

4000 

2292 

1.00 

2913 

MFO/LUNG 

30/48 

2000 

2238 

0.53 

1457 

41/48 

4000 

2369 

1.00 

2913 

MFO/LIVER 

16/48 

2000 

5155 

0.53 

1457 

40/48 

4000 

5475 

1.00 

2913 

DCM/LUNG 

30/48 

2000 

477 

1.09 

1457 

41/48 

4000 

1040 

1.00 

2913 

DCM/LIVER 

16/48 

2000 

426 

1.14 

1457 

40/48 

4000 

976 

1.00 

2913 

PATHWAY/ 

EFFECT 

NTP 

BIOASSAY 

RESPONSE 

MOUSE 

EXPOSURE 

DOSE 

GLOBAL83 

HUMAN 

INTERNAL 

GLOBAL83 

95%  LCL  at 
1.00E-08 

ppm 

mg/kg/day 

mg/kg/day 

GST/LUNG 

30/48 

2000 

1336 

1 . 2753E-05 

41/48 

4000 

2913 

GST/LIVER 

16/48 

2000 

1273 

3 . 2764E-05 

40/48 

4000 

2913 

MFO/LUNG 

30/48 

2000 

2751 

1 . 9424E-05 

41/48 

4000 

2913 

MFO/LIVER 

16/48 

2000 

2743 

3 . 1728E-05 

40/48 

4000 

2913 

DCM/LUNG 

30/48 

2000 

1336 

1.2753E-05 

41/48 

4000 

2913 

DCM/LIVER 

16/48 

2000 

1272 

3.2692E-05 

40/48 

4000 

2913 

The  modification  factor  is: 


2000  ppm  *  dose  surrogate  at  4000  ppm 


4000  ppm  *  dose  surrogate  at  2000  ppm 
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PATHWAY/ 

EFFECT 

HUMAN 

EXPOSURE 

HUMAN 

DOSE 

SURROGATE 

CST/LUNC 

GST/LIVER 

PP« 

2.88E-04 

4000 

2.88E-04 

4000 

ag/ltr/day 
1 . 7606E-06 
55.45 

3 . 8484E-06 
297.48 

MFO/LUNG 

MFO/LIVER 

2.88E-04 

4000 

2.88E-04 

4000 

2. 1215E-06 
1.21 

8  3471E-04 
1195.38 

DCM/LUNG 

DCM/LIVER 

2.88E-04 

4000 

2 . 8BE-04 
4000 

3.3210E  -05 
1045.96 

1 . 3067E-05 
1010.05 

PATHWAY/ 

EFFECT 

HUMAN 

EXPOSURE 

UNIT  RISK 
95X  LCL/ 
Target  Dose 

GST/LUNG 

GST/LIVER 

pp» 

2.88E-04 
2 . 88E-04 

9.89E-08 

1.57E-08 

MFO/LUHG 

MFO/LIVER 

2.88E-04 

2.88E-04 

3.56E-06 

8.74E-07 

DCM/LUNG 

DCM/LIVER 

2.88E-04 

2.88E-04 

9.89E-08 

1.57E-08 

DOSE  SUR. 
MOD 

FACTOR 

HUMAN 

INTERNAL 

TARGET 

DOSE 

Int/Mod 

ag/kg/day 

ag/kg/day 

2.26 

2 . 8571E-04 

1 . 2615E-04 

5.56 

2. 8571E-04 

5 . 1401E-05 

0.04 

2.8571E-04 

6 . 9598E-03 

0.10 

2 . 8571E-04 

2 . 7744E-03 

2.26 

2 . 8571E-04 

1.2615E-04 

5.56 

2.8571E-04 

5 . 1402E-05 

COMBINED 

UNIT 

RISK 

lung+ liver 

1.15E-07 

BODY 

WEIGHT 

1.45E-06 

BODY 

SURFACE 

4.46E-06 

BODY 

WEIGHT 

5.64E-05 

BODY 

SURFACE 

1.15E-07 

BODY 

WEIGHT 

1.45E-06 

BODY 

SURFACE 
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PBPK  method  -  Body  weight  (and  Body  surface)  -  PBPK  method 


combinations  (Andersen  and  Clewell,  1988): 


PATHWAY/ 

EFFECT 

NTP 

BIOASSAY 

RESPONSE 

MOOSE 

EXP. 

DOSE 

MOUSE 

DOSE 

SURROGATE 

DOSE 

SURROGATE 

MODIFIED 

CLOBAL83 

DOSE 

SURROGATE 

GLOBAL83 
95%  LCL  at 
1.00E-08 

ppm 

mg/ltr/dy 

mg/ltr/dy 

mg/ltr/dy 

mg/ltr/day 

GST/LUNG 

30/48 

2000 

314 

53 

53 

5 . 0456E-07 

41/48 

4000 

684 

115 

115 

GST/LIVER 

16/48 

2000 

1001 

168 

168 

4 . 3126E-06 

40/48 

4000 

2292 

385 

385 

MFO/LUNG 

30/48 

2000 

2238 

376 

376 

2.6538E-06 

41/48 

4000 

2369 

398 

398 

MFO/LIVER 

16/48 

2000 

5155 

866 

866 

1 . 0027E-05 

40/48 

4000 

5475 

920 

920 

DCM/LUNG 

30/48 

2000 

477 

80 

80 

7.6501E-07 

41/48 

4000 

1040 

175 

175 

DCM/LIVER 

16/48 

2000 

426 

72 

72 

1 . 8686E-06 

40/48 

4000 

976 

164 

164 

PATHWAY/ 

EFFECT 

HUMAN 

EXPOSURE 

HUMAN 

DOSE 

SURROGATE 

UNIT  RISK 
95%  LCL/ 
Target  Dose 

COMBINED 

UNIT  RISK 
lung+ liver 

GST/LUNG 

GST/LIVER 

ppm 

2.88E-04 

2.88E-04 

mg/ltr/day 

1 . 7606E-06 

3. 8484E-06 

3.49E-08 

8.92E-09 

4.38E-08 

5.55E-07 

BODY 

BODY 

WEIGHT 

SURFACE 

MFO/LUNG 

MFO/LIVER 

2.B8E-04 

2.88E-04 

2. 1215E-06 
8.3471E-04 

7.99E-09 

8.32E-07 

8.40E-07 

1.06E-05 

BODY 

BODY 

WEICHT 

SURFACE 

DCM/LUNG 

DCM/LIVER 

2.88E-04 

2.88E-04 

3. 3210E-05 

1 . 3067E-05 

4.34E-07 

6.99E-08 

5.04E-07 

6.38E-06 

BODY 

BODY 

WEIGHT 

SURFACE 
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